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Analytical mechanics of a spinless massive particle is generalized to that in which quantum phenomena 
appear. This generalization is accomplished with two keys. The first is the extended action of the 
diffeomorphism group on the covariant tensor fields of order two, with which the classical Hamilton- 
Jacobi eq. of the particle is modified. The other is the second-order Lagrangian L(q,q,q), with which 
the Lagrangian formalism consistent with the modified Hamilton-Jacobi eq. is constructed. In the thus 
obtained formulation, i) the particle picture is held at all times, ii) the discrete energy levels of the 
system of a particle in a confining potential are the same as those of quantum mechanics, iii) the 
particle distribution in an ensemble differs sHghtly from that of quantum mechanics, which however 
does not mean the incompetence of the present formulation because the precision of the experimental 
data up to date on the particle distribution is not enough to exclude the present formulation. 
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1. Introduction 

Despite the great success, serious conceptual difficulties remain unresolved in Quantum 
Mechanics (QM) (e.g., Ref. T chap. 17). One of them is that the particle picture cannot be 
maintained in (the Copenhagen interpretation of) QM; the so-called double-slits problem 
illustrates this difficulty. There have been attempts to recover the particle picture in QM. 
Notable is the Bohm's interpretation (B.I.), which allows a realistic particle picture in non- 
relativistic QM. However, there arise other difficulties in B.I. Most serious would be that it 
does not work for a relativistic spinless particle (for B.I., see Ref. 2 and references therein)0 
In Ref.jSi Faraggi and Matone made a breakthrough along the line of B.I. They derived 
the field eq. used in B.I., which we call the Quantum Stationary Hamilton-Jacobi Equa- 
tion (QSHJE), modifying the classical Hamilton-Jacobi equation (H-J eq.) with what they 
call 'the Equivalence Principle', which states that all physical systems can be connected by 
a coordinate transformation to the free one with vanishing energy. (In classical analytical 
mechanics (CM), if the energy E of the system vanish: E — Q, then §o(->c) — const., where 
So('*^) is the reduced action. Accordingly, we cannot transform the system to the £ 7^ 



" Although talking about the relativistic trajectory goes beyond the scope of the present paper, we can show that 
the argument against the possibility of introducing the trajectory to a relativistic spinless massive particle in B.I. 
does not apply to the present formulation; see 
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system with the time independent coordinate transformation.) To be noted are i) they ob- 
tained the QSHJE not from the Schrodinger eq. (SE) of QM (in B.I., it was obtained from 
SB) but from the H-J eq. of CM, and ii) discrete energy levels appear without recourse to 
the axiomatic interpretation of the wave function of QM. They appear, in one-dimensional 
(1-D) space, from the requirement that e^'SoW/fi^ where So(x) is a solution of the QSHJE, 
be continuous at x = {±°°} in the projective geometrical sense. Their formulation however 
lacks, as a CM, the Lagrangian formalism in which the trajectory is determined. 

In Ref. 4 Bouda determined the 1-D Lagrangian L consistent with 1-D QSHJE in the 
form: L{q,q,q,q) and derived the 1-D Euler-Lagrange eq. (E-L eq.) from it. Bouda de- 
termined them considering that the L should depend not only on q, q but also on q,q, - ■ ■ 
(q — dq/dt and so on) because 1-D QSHJE is a third order ordinary differential eq. 

Hinted at by their works, we formulate a quantum theory of a non-relativistic spinless 
massive particle as a generalized, or modified, CM. It does not suffer from the difficulties 
of B.I. (gli. 

We outline our formulation. 

Our formulation consists of the H-J formalism and the Lagrangian formalism (the latter 
is transformable to the Hamiltonian formahsm; [Appendix C| l. The H-J formalism is based 
on the QSHJE: 

where 7? e M and §o G are unknown functions. As a partial differential eq., the solution 
gives a momentum {V$^field. In the region where —h^V^R/lmR, which is called the 
quantum potential (QP),'^ is not ignorable, the VSq can be wave-like if we give suitable 
initial conditions. The Lagrangian formalism is based on the semiclassical Lagrangian (in 
principle, we use, as coordinate variables, x and its variants in the H-J formalism, q and its 
variants in the Lagrangian formalism): 

2 2m 

where e(q) << 1 is determined, in principle, to make §o(x) — Et — /" 'Lt/f be true for a 
solution So(x) of ( II. lb . Accordingly the i(q) makes it possible to determine the trajectory 
even in the interference region. The q^/q^ makes the classical trajectory wavy. The second 
term of the RHS corresponds to the QP appeared in ( II. 11 1. The Lagrangian formalism is 
effective at classical and semiclassical levels; at the quantum level, our formulation has no 
Lagrangian ( m.2\ . 

We derived (II. Il l modifying the classical H-J eq. with the assumption that the action of 
the diffeomorphism group on the covariant tensor fields of order two is extended. We note 
i) although (II. lb is obtained from the SE of QM replacing the wave function ^ in it with 
^giSo/fi taking the real and imaginary part of the resulting eq. (e.g., Ref.|2|l, we do not 
use this method since, for us, QM is not a starting point but a goal, ii) as was mentioned 
above, Faraggi and Matone derived (II. lb with 'the Equivalence Principle'. 

We determined our Lagrangian to the form of (11.2b as a quantity which makes §o(x) — 
Et — p ' Ldt true for a solution §o of ( II. lb . We note i) the 1-D Lagrangian of Bouda SI 



\me- ''^^'''" ^-Viq), (1.2) 
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differs from the 1-D form of (11.21 ): there are some differences between his and our method 
used to derive the Lagrangian, ii) the Lagrangian formaUsm is constructed on the higher- 
order-analytical mechanics, for which see Ref.|5]|6]and references therein. 
Quantum phenomena appear as follows. 

• Uncertainty relation: In a 1-D square potential well, the momentum p := VSq constructed 
from the solution of the 1-D version of ( 11.11 ) undulates. Moreover, the narrower the well, 
the larger the amplitude of the undulation, which corresponds to the uncertainty relation 
between position and momentum of QM. 

• Energy and angular momentum quantizations: In 1-D space, energy quantization appears 
from the requirement that g^'SoW/fi jg continuous at x = {±°°} for a solution 2>q{x) of (II. 11 1, 
which is equivalent, if the V is the confining one, for the second order hnear differential eq. 
corresponding to ( fTTT ): -{h^ /2m){d^^' / dx^) + {E~ V)^ = to have a L^-solution. Since 
this eq. is the same as SE, we see that the energy of the system is quantized to the same 
levels as those of QM.Elln 3-D space, the energy, the angular momentum too, is quantized 
by the 3-D version of the 1-D requirement. 

• Tunneling: In our formulation, the kinetic energy of the particle depends not only on q 
but also on q and q. Accordingly, even if the potential hill is higher than the kinetic energy 
of the incoming particle, the velocity can be non-zero real: q^ > at the top of the hill. 

• Interference: We consider the double-slits experiment for concreteness. In which, the in- 
terference pattern on the screen appears as a result of undulating trajectories. The trajectory 
in the interference region is determined using the formula connecting VSq and (q,q,q): 

V§o(x) — mq+ (/!^/m)(q/q^) H |q=x (see ( 12.36b ). We regard this formula as a trajectory 

determining eq. after inserting the VSo constructed from the solution of ( II. lb into the LHS. 
The effect of the double-slits manifests in the functional form of the VSq. We cannot use 
the E-L eq. for the present purpose because we cannot give the exact form of e(q) in ( 11.21 ). 

We note above reasoning for 1-D energy quantization is essentially the same as that in 
Ref.[3] We however repeat it for completeness. Other reasoning would be found nowhere. 

Our formulation predicts that the density distribution (dens, distr.) of massive particles 
in an ensemble deviates from j'Pp (by a few percent in the experimental setting of '> \2>3\ , 
which however does not contradict existing data supporting QM because the precision of 
the data is not enough to exclude our prediction. For example, the exper imental data of the 
dens, distr. of electrons Fresnel-diffracted by the electron biprism,l21SJ9|yYjjj(;jj supposed 
to support QM's prediction, does not exclude our formulation because the experimental 
error (here we mean the difference between the measured and quantum mechanically cal- 
culated dens, distr.) by more than several percent ^ is larger than the expected difference 
between QM's and our prediction by a few percent. Thus, we see that our formulation has a 
chance to be verified by experiments. (Casual inspection makes us feel that the data support 
our formulation; see the last paragraph of ^^3.31 ) 

The organization of the present paper is as follows. In ^ ^2.11 we derive the QSHJE. 
In ^2.2\ we construct the Lagrangian formalism. In f|3] we show that above mentioned 
quantum phenomena appear. In particular, in ii3.3l we show that the dens, distr of diffracted 
electrons differs from that of QM. In [0] we describe how our formulation differs from B.I. 
In appendix, we put background materials and lengthy formula derivations. 
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2. Formalism 

2.1. Quantum Hamilton-Jacobi equation 

We first modify the relativistic H-J eq. assuming the extended diffeomorphism action. Then, 
we take the non-relativistic Hmit of the modified H-J eq. 

The extended action of the diffeomorphism group of 1-D projective space MP': 
Diff(RP') upon the covariant tensor fields of order two on RP': 3^2 (KP') is possible; see 
Ref. 10 and |AppendixAl The (not-extended) action T of Diff(RP^) on J2(RIP') is given as 



T : Diff(Mpi) X J2(KP') ^ :f2(RP'), {y'\<P<j)) ^ Ty-i(j){y) = (j){y{x)){dy/dx)^. While 
the extended action (we write it as T'frf,.-! y-i)) is given as 

Tidy-^,y-^)Hy) = (^)20(3'(x)) + a{y(x),4 , (2.1) 

where a = const. G R and {3;(x),x} = y — j (^)^ is the Schwarzian derivative (/ = dy/dx 
and so on). We call this action the extended action of the diffeomorphism group or the 
action of the extended diffeomorphism group. 

The extended action of the diffeomorphism group of 4-D Minkowski spacetime on 
the covariant tensor fields of order two on M"*: Tf^pCMf^) is also possible (the 3^2 (KP') is 
identified as a restriction of Tf^j^ (Jvf*) to 1-D space). Indeed, we constructed the extended 



action 'on' the H-J eq.; see Appendix B We consider the diffeomorphism of M as the one 
between patches not including infinity. We therefore do not care about the global topol- 
ogy of the base manifold. With the extended action, the H-J eq. (our Minkowski metric is 

diag(l, -1,-1,-1)): 



where § is the relativistic action, e is the electric charge of the particle, m is the mass of the 
particle, c is the velocity of light, and (A*', A) is the electromagnetic potential, is modified 
to what we call the relativistic quantum H-J eq.: 

(5|.S + %)(<?''§ + -A") = m2(c2 + e), Q=tM^ (2.3a) 
c c Rm'^ 

^m/=0, f^^^idf^S + '-A^), (2.3b) 
m c 

where fi runs over ~' 3, and /? is a function determined solving the equations. 

Taking non-relativistic limit, that is, setting S = S,„- — mc^t (S„,. is the non-relativistic 
action) and ignoring • /c^ terms (we consider only the case: A' = A^ = A^ = and is 
time-independent; in this case, R is time-independent and S„, is decomposed as S„,-(x,r) — 
§o(x) — Et, where So is the reduced action), we have the non-relativistic quantum H-J eq. 

where ; runs over 1 , 2, 3. If space dimension is one, ( 12.41 ) is simplified, using the formula in 
footnote |x] to 

2L(5)' + ^{^<"^>+^W-^-0, (2.5) 
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where Vix) = eA^{x). We call and (US 3-D QSHJE and 1-D QSHJE respectively after 
Ref. 300 We call h^V^R/{2mR) and h^{§,Q,x)/{Am) 3-D and 1-D quantum potential (QP) 
respectively. Replacing q" — > x, ^ §o, q'^ — > e^'®"/^' in the cocycle condition: {<?' ,<?"} = 
{dq'' /dq"f{q'\q''} + {?^?"}, which is equivalent to d{g-^h-^) = dg-^ + T^-idh'^ on 
pl30] (these replacements are allowed because x M- §o{x). So g^'So/s^ gj^. jg^gj 
local diffeomorphisms; see footnote|v)l, we rewrite i2.5\ as 

V(x)-E = ~—{eT^^o,x} . (2.6) 

Using the theorem: 'Let the two linearly independent solutions of d^z/dx^ + u{x)z — 
be zi,Z2- Then, {zi/zi^x} = {{Azi +Bz2)l{Cz\ +Dz2),x} = 2u{x), where A, - ,D = 
const, e C and AD -BC ^ 0.' (e.g., Ref.|10j|), we solve the 1-D QSHJE as 

where W and 'P^ are two linearly independent solutions of the corresponding (having the 
same V{x) and E) second-order linear differential eq.: d^z/dx^ — (2m/h^){V (x) — E)z = 0. 
Since this eq. is equivalent to the SE of QM, we henceforth call this eq. SE. Keep in mind 
however that in our formulation SE is used only as a mathematical tool to solve ( I2.5l l (or its 
3-D version ( I2.4l l as is seen on p|7]i. We do not give any physical meaning such as 
represents the probability density of finding a particle at x ' to or M^^. Physics is discussed 
with §0- 

The momentum p : — dSo/dx is obtained differentiating both sides of ( I2.7l i with respect 
to q as 

„ i HAD-BC)W -^^^^^ 

where W is the Wronskian: W ~ _ const. ^ 0. We require p to be real 

because p is a classically observable quantity. We make ( 12.81) a little readable. We set, 
without losing generality, 'P, S M and AD — BC = ±2/ for convenience. Then, with /j = 
i{AD - BC) /2AC = Tl /AC g M\{0} and Z2 = (AD + BC) /2AC G M, (ESl is transformed 
as follows with p being kept real: 

^ dx (^0 + /2^)2 + Z2>I/2 • ^ ■ ^ 



In accordance with the appearance in (12. 6t of the Schwarzian derivative, which lives in 
1-D projective space MP', we put a requirement on the So that, as a physically acceptable 
solution, 

the e^'^o/s is continuous on R := RU {±°o}(« RP^). (2.10) 



As was mentioned in f[T] Faraggi and Matone derived \2.5\ in Ref. 3 through a different path from ours, through 
what they call 'the Equivalence Principle'. 

■^In writing this and following paragraphs to the end of this subsection, con'esponding pai1s of Ref. 3 were helpful. 
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Since e^'^o/'' is always continuous (even differentiable) on R as is seen from ( |2.9l l. this 
requirement is essentially the continuity 00 requirement at x = {±00}. Since the e^'So/s 
is continuous at x = {±00} if and only if lim;r^+„oSo(x) — lim;c_j._„c§o(-^) — nh/2 {n = 
±1 , ±2, • • • ), (12.10b is equivalent to the requirement: 

f+°° nh 
d§o~ / pdx ^ — , (2.11) 

which corresponds to the Bohr-Sommerfeld quantization condition: f pdx — f d§o — nh. 
We consider that the e^'So//; continuous at x = {±00} for n = ±00. If we use ( 12.1 11 1 instead 
of ( 12.101 ). we need not refer to the neighborhood of infinity; if the physical system under 
consideration is of atomic scale, integration interval of ( 12.1 11 1 is taken to be much less 
than, say, 1/im. Thus, although our formulation is constructed on ]RP\ the formulation is 
indistinguishable from the one on M on which CM and QM are considered to reside. 
In the case of a confining potential, ( 12.10b is equivalent also to the requirement: 

9 

The corresponding SE has a L -solution, (2.12) 

which is seen as follows. Inserting ( 12. 9b into ( 12.1 lb and using a formula: J[{f' {x)g{x) — 
f{x)g'{x) )/{ f^{x) +g^(x) )]dx — arctan[/(x)/^(x)], we rewrite ( 12.1 lb as 

V^x^r^, , ,.x = -;.arctan^:^r^^, (2.13) 

which is equivalent to require for the ^^/^ to be continuous at x = {±°°}- In the case of 
a confining potential, the ^^/'P is continuous at x = {±°°} if and only if corresponding 
SE has a L^-solution, which is seen from i) ^ iii) below. Let ^'c and be two linearly 
independent solutions of the 1-D SE having a confining potential. Then, 

i) The ^'c can be an L^-solution only for discrete values of E. For such £"s, the other 
solution diverges forx = ±00.^^ 

ii) If the is a L^-solution, ^^/^^ is continuous at x = {±00}, which is seen from 



l'^(x) =c^(x) / '^-^-{q)dq + d^'{x) , (2.14) 

■to 



where c,ii G C (c 7^ 0) are consts. (the formula ( 12.14b is obtained solving W = ^>'^>'^ — 
^ ~ const . regarded as a differential eq. for with the variation of constants method), 
iii) If none of I'c and is an L^-solution, ^'^/^c cannot be continuous at x = {±°°} 
because and 'Pf behave forx — > ±00 as follows (Ref.[T4lchap.3). 



A function f{x) e R on R is said to be continuous at ,ro £ R if 

it-vo [-lim_j^j^/(.r) for lim.v^.v„ /(x) = +~ or - 00 . 

■^We cannot require differentiability at x = {±00} as is seen from the following example. Consider the particle 
moving in a non-confining potential. Let the energy of the system be E and the potential approaches to a con- 
stant value V- <E for a: ^ -°o and y+(^ V-) < E for x +°o. Then, Um.,.^_o„3So/5-v = \/2m(£-y_) ^ 
y'2m{E-V+) = limi-^+„ 3So/ dx . 
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For X ±00 and 0; for x +00, and ^ ±00. Or, 

for x-i' -oo,Wc-)'0 and ^ ±00; for x +00, ^ ±00 and ^ 0. 
Thus, in the case of a confining potential, ( 12.12b is equivalent to (12.1 11 1 or ( 12.10b . 

If space dimension is three, the solution set (R, So), where R,§o E M, of 3-D QSHJE 
(12.4b is constructed from the solution \j/ of the corresponding SE as = Re'^'^l^ since ( 12.41 ) 
is equivalent to the SE with = Re'^'^l^^ (e.g., Ref. (2). The momentum p := VSq is given 
as 

.^^^V^^-^ (2.15) 

where 1/ is the complex conjugate of y/. Of course, (12.15b is valid also in 1-D space. In 1-D 
space, the p obtained from ( 12.151 ) setting 1/ = ci^ + C2^^ is equal to the p obtained from 
(12.8b with A = ci,B = C2,C = ci and D = c^. The requirement imposed on 3-D So is given 
in El 

Example We consider a 1-D free particle of energy £ > 0. 

If £ > 0, we take ^ = coskx and = sin^x, where k — y/2mE/Ti, as two linearly 
independent solutions of the corresponding SE. The momentum p is given, from (12. 9b . as 

^ ::m 

(sin kx + h cos kx)^ + 1\ cos^kx 



If /i —^l and I2 = 0, we have the classical value: p = ±\^2mE. 

If £■ = 0, we take ^ = 1 and — x, where 1 is a constant having the same dimension 
as X. The momentum p is given as 

_ ^So _ -hhl 

dx (x + /2l)2 + (/il)2 • ^ ' 

The e^'So/fi jg continuous at x = {±00} since we have lim;r^+ocSo(x) — lim,v^_ooSo(x) ~ 
±hl2 from (ITTtI i. The So constructed from e^'So/fi ^ [Ax + B\) / {Cx + D\) is therefore 
physically acceptable with A ~ D which make p real. 
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2.2. Lagrangian formalism 

From the solution So(x) of the QSHJE, we learn the value of the particle's momentum 
p := V§o at X. From the QSHJE, however, we cannot see whether the particle comes to x or 
not, nor see we how p is related to (q, q, • • • ) of the particle (VSo(x) — mq|q=x is no longer 
valid). To see these, we construct the Lagrangian formalism. 

We consider, based on a philosophy of CM, the eq. of motion exists as an ordinary dif- 
ferential eq. (ODE) even at the quantum level. Although not all ODEs are obtained with the 
variational method from the Lagrangian El we search for the eq. assuming that our sys- 
tem has the Lagrangian. If not, we will see contradictions somewhere in the eq. searching 
process. Indeed, we will find that the Lagrangian does not exist beyond the semiclassical 
region, at least within our theoretical framework. 

Our semiclassical Lagrangian cannot be a function of q and q as in CM. We show the 
reason. In the Lagrangian formalism, the p is a function of the trajectory q(f). That is, the p 
at the particle position q(f) is given as p|q(r) = • • • )l' " ' (*!)■ Let us consider a 1- 

D free particle. Assume (to get a contradiction) that our Lagrangian of which is a function 
of q and q: L{q,q). Then, the E-L eq. is of second order. The trajectory is therefore fixed if 
two integration parameters, say, particle position qo at time Iq: ^o(fo) and the energy E, are 
given. The momentum p at any point on the trajectory is also fixed from (* 1 ). However, the 
formula ( 12.16b tells us that, even if qo and E are given, the p at, say, is not fixed, which 
implies that the trajectory is not fixed with ^o(fo) and E only. 

If the E-L eq. is of fourth order, the above problem is resolved because four parameters 
(x,E,li, I2) necessary to fix the p of ( 12.161 1 correspond to four integration parameters of the 
E-L eq. Thus, we introduce the higher-order Lagrangian L{q,q,q, ■ ■ ■), from which fourth 
order E-L eq. is derived. Or we assume that the system described by the QSHJE has the 
Lagrangian with §o{x)—Et = J^'' L{q,q,q, ■ ■ ■ )dt ■ ■ ■ (*2), where §o(x) is a solution of the 
QSHJE. 

We note on a seeming inconsistency in the (*2). In the second (similar for higher) order 
Lagrangian formalism, the reduced action So(x,x) is obtained as a function of x and x from 
§o(x,x) ~ Et ~ f^'^'' L{q,q,q)di ■ ■ ■ (*3), where x is the velocity at the end point x (see 
(IC.ll l. see Ref. 22 §9). Whereas the reduced action §o(x) obtained from the QSHJE is a 
function of x. This seeming inconsistency is understood as follows. In (*2), the integration 
is carried out along the trajectory. Since the trajectory, which is determined by the E-L eq., 
depends on initial (ini.) conditions, the So(x) too depends on the ini. conditions of the E-L 
eq. — we may write §o{x) as So(x; ici ,ic2, - ■ ■ ), where ici , • • • are ini. conditions. If we take 
x(f ) and x(f ) as parts of the ini. conditions, the (*2) is written in the form of (*3). Thus, 
we see that (*2) and (*3) are essentially the same; In (*3), parts of the ini. conditions x(r) 
and x(f ) are explicit, whereas in (*2), ini. conditions other than the x(f ) are implicit. 

In the following, we determine the semiclassical Lagrangian of our formulation making 
use of the higher-order CM ( [Appendix C| l. First, we determine the Lagrangian of a free 
particle in 1-D space. Then, we generalize it to that of a particle in a 3-D potential field. 

We expand the would-be 1-D free particle Lagrangian L — L{q,q,q,q,- ■ ■) in power 
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series of h: 



1 



We assume that the physical parameter (the parameter having dimensions) appearing in 
each term /i ,/2, • • • is only mass m. Then, from the dimensional analysis, we see that /j.'s 
are restricted to 



/i 



ciq 
4 

where ci,C2a,- 



/2 



C2a q 



or/2 



C2b q 



fk 



Ck 



ot" ' q 



3k-2 



m q-' m q^ 

are dimensionless parameters independent from q, ■ 



: 3 4 • 

^ , -r, 



(2.18) 



2 q 



2 y:2 



q , C2,fi'- q C2bh q 



q' 



q^ 



I 

k=3 



ci,h q' 



, q. Thus, we have 
(2.19) 



The second and higher terms of the RHS as a whole is the Lagrangian side counterpart of 
the QP appeared in the QSHJE ( I2.5l l. We call it too the QP. The ci,C2a, • are determined 

by /], I2 and k = \/2mE/h of ( 12.161 1. Although ci ,C3,C5, terms are not invariant under 

time reversal, we keep them intact for the moment. The ( |2.19t is not defined for q = due 
to the q in denominators, which however causes no problem because the ^ of a solution of 
the E-L eq. never takes value zero unless ^ = (Appendix El. 

The E-L eq., J-O momenta /^p), P{3), the energy E obtained from ( 12.191 ) are 
( Appendix C| and [Appendix^ 1 



— mq- 



2{3c2a+C2h)h^ 



q 8qq lOq' 



{k-2)kcf-^q^- 



q^ 

-kif^^'q 



2(3k^ 



^3k-2 



q° 

~ 2)kq'' 

^3k^l 



^ {k^l)ckh'' 

fc=3 

^q {3k-2){3k-l)q' 



■k+l. 



(2.20) 



P{2) 
P(3) 



dL d dL (f dL 
dq dt dq dfi dq 



k=3' 



2{3c2a + C2b)Ji q , 4(3c2a + C2b)Ji q 



2 ;;2 



m q^ 
y,k-2y, Uy,k 



k{k-l)ct-^q {3k-2){k-l)ct 



q 



3k-2 



dL d dL cifi ^{3c2a + 2c2b)h^ q 
dq dt dq q m 

dL _ C2cfi^ 
dq mq^ 



kc^ti q'^ 



k -k-l 



q^ ' q^^ ^ 



q-' 



(2.21a) 
(2.21b) 
(2.21c) 



2 ^2 



E = P(\)q + P{2)q + P(3)q- ijmq 



2{3c2a + C2b)fl q , 5{3c2a + C2b)fl q 



Ckfi 



k=3' 



vfk 1 



m q^ 

k-2y, (IK 1\y,k^ 



k{k -l)q''-^q {3k-l){k-l )q^ 



q 



3k-3 



q 



3k-2 



(2.22) 
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Note i) the E-L eq. is fixed to the fourth order only from the dimensional analysis, ii) the 
c\ does not appear in the E-L eq., p^y^ and E. 

The momentum pel d§(){x)/dx obtained from (12. 19b is 

PEL^Pil)+Pi2)^+Pii)^ 

q C2ji^ ii , {2>C2u+2c2b)h^ if' , ^ kckh'' (f 

= P(i)+c\h— + — + — + Lt:T^3FT' "^2.23) 

^ ' q^ m q* m q^ m" ^ q'^" ' 



which is obtained as follows. We know (Appendix C i: 

§ = J L{q,q,q,ii)dt , E = p(Y)q + p(^2)q + P{i)ii - L . (2.24) 
From ( 12.241 1. along the solution path, we have 

d$ = Ldt = {p(i)q + P(2)q + P{i)q -E)dt 

= P{i)dq + P{2)dq + P{3)dq - Edt = {pi^i-^ + P(^2)^ + P(j-^^)dq - Edt . (2.25) 

q q 

Comparing ( 12.251 ) with d§ = c/Sq — Edt = pdq — Edt of CM, we see that pel = P(\) + 
P{2)q/q + P{3)q/q- Thus, we have ( 12.231 1 (the p^ij, which is const, from ( IC.6b . is left intact 
for later convenience). 

For the E-L eq. ( 12.201 1 to be the eq. of motion, it has to be 3c2„ + C2b < 0. Otherwise, 
a wavy linear trajectory required for the pel to correspond to the p of (12.16b cannot be a 
solution of (12.20b . We call p of ( 12.16b Pqhj henceforth. 

If C2 := 3c2a + C2b < 0, for any semiclassical initial conditions, that is, for {q{to),- ■ ■ ,q{to)) 
satisfying ^(fo) ^0, |^(fo)| « tn\a(to)\y {hy/ -2c2) and |^(fo)| « m^\q{tQ)\y {~2h^C2), 
the solution of (12.20b is written as|3 

q{t) = vf + ^ {a„cosn(Ot + b„ sinnof) = v'f + ^ A„cos{n(Ot + 9„) =: vt + f , (2.26) 

where |v| — const. > and |//v| << 1; an additive const, was set to zero. The oscillat- 
ing part is written as a Fourier series ( [Appendix E| i. The infinitesimal (inf.) order of A„ is 
(f/v)" ( [Appendix E| i. The v, ai, bi (Ai, 0i) and the omitted additive const, work as initial 
conditions of the fourth order E-L eq. The (O is determined to (O = mv^/ {h^/—2c2 ) from 
the E-L eq. truncated at (//v)'. The (A2, 62), (A3, 63), • • • are determined recursively from 
(Ai , 01 ) and v. For example. 



A2Cos(2wf + 02) = -1/1 % — At cos(2a)f + 201 +^) , (2.27) 

8c| V 



^ We warn the reader that the w in 12.26) is not the (0 appearing in QM as the angular frequency of a massive 
particle, say, in a formula: E = p^/2m = ha). Let the o representing the angular frequency of the massive particle 
be (Oq, the £0 in i2.26\ be (O^. In non-relativistic QM, from p = mv = hk and t' = dwQ/dk, we have E = /2m = 
hcOQ. While in our formulation, approximately, although p = mv = tik is common, v is expressed as t' = ffli/i. We 
therefore have E = p^/2m = Tia>ij2 (see )F. 14) ). 
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2(-2c2)^ V 6 V 

+ (4 + ^)(— cos(3«f + 30:) Z£^(^)2ai sin(3a)f + 30i) , 

24 32c| V 12(-2c2)2 ^ 



where a> ~ mv^ j (h\J —2c2 ); the ^ is the angle satisfying cos = ( 2/^/2 /\/"*" ~ '^3/(^''2) 



and sin^ = -c^/(8c^) . 

If the initial conditions are not of semiclassical, the solution is not necessarily wavy 
linear. Assume that initial conditions {q{to), ■ ■ ■ ,^(fo)) satisfy q{to) > and P{i) — Vipfi) ^ 
2mE) < 0. (Note, if initial conditions is semiclassical, — V {pfij ~ 2mE) ~ pji) > and 
{p{i) — %/ (p^i) ~ 2mE))/m is the minimum velocity ( Appendix E| l.) Then, for t — > +00, the 



q asymptotically approaches some const. > because q cannot cross ^ = (Appendix E 1. 

The Lagrangian ( 12.191 1 in which at least one of ci and C2a is not zero and the one in 
which ci — C2a = describe different systems in the sense that the relation between k and 
(O is different: O) = 2k\' for the former, O) = kv for the latter, which we see comparing pqhj 
and Pel at the lowest inf. order of oscillation after inserting (12.26b into x of ( 12.16b and q of 
(12.23b . In the present paper, we study the latter because we consider the latter is likely to 
describe the nature for reasons below, i) If ci 7^ (in general, if c„ 7^ 0, where n is odd) the 
Lagrangian is not invariant under time reversal, ii) If we require our free 3-D Lagrangian 
to be rotationally invariant as in CM, it is natural to consider that ci = C2a = (in general, 
c\ = C2u = C3 = C5 = • • • = 0) because the c\- and C2a-terms (in general, cj-, C2a-, C3-, C5-, 

terms) of ( 12. 19b has no natural such 3-D form, which is contrasted to other terms having 

such 3-D forms with q^ ^ i\ - ii~ q\ + q2+ q\ , etc. iii) If C2a 7^ 0, the Lagrangian is of the 
third order. However, because the fourth order E-L eq. is obtained from the simpler second 
order Lagrangian, we have no reason to adopt the third order Lagrangian. 

We determine C2b {c\ = C2a = 0) comparing Pqhj and pel- For the present purpose, we 
write Pqhj or (12.16b in the form: 

^ (2.28) 

(sin^x + e2Cosfcx)2 + (l+£i)2cos2fex ' 

In accordance with the present semiclassical setting, we replaced /j — ^ — (1 +£1) and 
h £2, where £1,62 << 1. Inserting q{t) of ( 12.26b into q of pel ( 12.23b and x of pqhj 
(12.28b . and comparing pqhj and pel at the lowest inf. order of oscillation or at e\ where 

£ = W£^ + £2, we have C2b together with fli,^i, CO and v as the function of £1, £2, and £ 



( jAppendixF] ). The results are: co — kv, v — y/2mE{l - £ + 0(£^))/w, C2h = -(1 - 2£ + 
0(£^))/2, where v > is assumed. The oscillation amplitude (ai,bi) are given as func- 
tions of £1, £2, and E in ( IF.16I) . dF. 17b . The inf. order of f/v with respect to £ is 1/2: 



^ It seems that we can construct rotationally invariant 3-D forms of the C[,C2a,C3, terms of )2. 19) if we take 

the absolute value of them. This idea does not work because the resulting E-L eq. is unacceptable. For example, 
the E-L eq. obtained from L = mq^ /2 + C2ah~\q\/m\q\^ is mq^ {icialf' I m)(2q / cf — \6qij/q^ + 20if /q^) = 0, 
where ' — '-sign for {q > 0, q > 0), '+'-sign for (q > 0, ij < 0), which does not give a wavy linear trajectory as is 
seen from the truncated eq.: mq^6c2ah^q/mq^ = 0. 
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f jv G 0(e'/^). The 1-D free Lagrangian is written as 



1 2 ^ (1 -2e + 0(e2)) ^' 
—mq 

2 2/n q' 



mi - "V^ --y-vc;; ^_^ Q(^^3) ^ (2.29) 



The C2i is fixed not only by the (ei,e2,-£') but also by the {q{to) ,q{tQ) ,q{to)) or by the 
(v,ai,bi) of (|Z26] | since they are related through dRlJl l, (IRT6T i, (IrTtT i. and ( |Z26] |. The 
dependence of C2fo on {q{to),q{to),q{to)) or on (v,fli,/7i) is admissible (Appendix D i. 

No value of C3 makes pqhj = /"el at the inf. order e^/^ of pqhj ( Appendix Gi, from 
which we conclude that there is no Lagrangian beyond the semiclassical region, at least 
within our theoretical framework. 

From ( |2.29t . omitting O(e^) and noting e = const., we obtain 

(E-Leq.) = ^^^1^ (1 - ^ + + 0(;.^) (2.30a) 

m \q^ q-' q° / 

p „^q-^ Z(^_^) + 0(ei) , p(2) = ^ --,+0{e--) (2.30b) 

^ ' m q^ q* ' m cf 



1 .2 ^2(1^ 5^2 ^ 

-mq {ttti- — 

L m Lcf q^ 

q . ;i^(l-2e) 3^2 ii 



E = -mq^-^—^{^^-^) + Oie--) (2.30c) 



3 , 



PEL = Pil)+P(2)-^mq-^ i(J^_^) + 0(e2) , (2.30d) 

^ ' ^ 'q m q^ q* 

where we have set C3 hidden in 0(/j^) of ( 12.291 ) to zero considering the Lagrangian is 
invariant under time reversal. A benefit of which is that we have 0(e^/^) rather than 0(e' ) 
in ( l230bb and ( l230c] i. Note 0(e^V2) = (3[[f/vf). 

Next, we determine the semiclassical Lagrangian of a 3-D free particle corresponding 
to a 3-D momentum: 

PQHJiW (sinkx + £2,COSkx)2 + (l+£i,)2cOs2kx ' ^ ^ 

where i— 1,2,3. We assume, from ( |2.19l l, the 3-D Lagrangian takes the form: L — mq^/2 + 
c^f /i^q^ /mq'* + (/z"^ ) , where q = (^i(f),^2(f).^3(f))'q^=q-q = ^i^i+?2^2 + ^3^3,and 
q^ = (q • q)^; likewise for q^. We set c^^ corresponding to C3 to zero (cj^ = 0) because there 
is no natural 3-D form corresponding to the C3-term (footnote [g]l. Then, in a similar manner 
as in the 1-D case, the is determined as the function of (£i/,£2/,£',) in the form: 



3D . 



my 



2„4 _„,2/„2 I „2 , „2\2 



'^26 



-m\v\ + v^2 + ^'i) 



2h^C0^ 2{\'\\/2mE\ + \'2\/2mE2 + V3^/2mEo,Y 



1 1 +0(e'), 



where co = t!=i kun, h = y/2i^i/h, Vi ~ (V2^/m)(l - ^ Ef- + E^ ), £ = l^ti Vi^ 
£2,) ( | Appendix!? I. 
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The semiclassical Lagrangian of a 3-D particle in a slowly changing potential field V (q) 
is naturally considered to be 

,.1 q2_^!(l_lMH;_v(q) + 0(,3), (2.32) 
2 Zm if 

where e(q) << 1, which is left undetermined, is a slowly changing function. 

The Lagrangian ( 12.29b and ( 12.321 ) is regular, that is the Hessian does not vanish 
{d^L/dif- ^ for ( IZ29t . \d^L/dq'dqj\ ^ for ( |232| |). a consequence of which is that 
the corresponding E-L eq. is fourth order and determines a unique vector field on the third- 
order tangent bundle T^M. ^ 

The 3-D free Lagrangian, or ( 12.321 1 with e(q) = const, and V{q) = 0, is translationally 
and rotationally invariant. The J-O momentum p^jj is the Noether current asocciated to the 

invariance of the Lagrangian under spatial translation. ^ 

From (|2.32t , ignoring 0{h^) and e(q) << 1 (if interference is involved, we cannot 



ignore e(q); see the paragraph after next), we obtain (Appendix C i 



dV{q) .. Jf- / qi 4^,qq 8^,qq 4^,qq 2^,qq 



(E-Leq.) = m^, + — (-— .g . .g . .g . .g 

oqi m \ q* q" q" q" q" 

^ 12g,(qq)(qq) 24g,■(qq)^ 33) 
qS q^ / 

^dL d dL ^ tP- ,qi ^ Iq^qi 4qq^,- . (2 34) 

dqi dt dcji ^' m q"* q^ q^ m q"^ 

The 3-D momentum p = VSq is obtained as follows. From S = / Ldt and E = P(i)q + 
P(2)q — i-, we see, along the solution path, c/S = dSo — Edt = Ldt = (P(i)q + P(2)q^£^)'^f = 
iP{i)i+Pi2)i qi/<ii)dqi + {■ ■ ■)dq2 + {P(i)3 + P{2)3 q3/qi)dq3-Edt. We therefore have 

= P{i)i + P{2)i— ^ mqi + — [ — + — z —j T4— lq=x ■ (2.36) 

axi ^ ' qi in \q* q" q° / m q*qj 

Although variables of 3-D free QSHJE separate on the Cartesian coordinates, those of 
the 3-D free E-L eq., or ( |2.33l l with V{q) = 0, does not, which however does not mean 
that two formalisms are inconsistent. The criterion of the consistency is whether physical 
quantities defined in both give the same values in both. In the present case, those defined 
in both are £0 and d§o/dxi. They give the same values up to e'. The two therefore are 
consistent up to e ' . We cannot discuss the consistency at higher orders of e because the 
Lagrangian formalism is not effective there (pfT2]i. 

In the interference region, we cannot ignore e(q) of ( I2.32l i if we use the E-L eq. to de- 
termine the trajectory. However it is difficult to determine e(q). We overcome this problem 



The author does not know whether i-dnection component £", of E is defined in the Lagrangian formalism or not. 
At least, there is no formula to express Ej. 
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bypassing the E-L eq., or reinterpreting ( 12.361 1 as a trajectory determining eq. Originally, 
( 12.36b gives dSo/dxi from a given q(r). However, since the dSo/dxj on the LHS is given 
also as a solution of the QSHJE, we regard ( 12.361 ) as a differential eq. for q(r) after insert- 
ing a solution of the QSHJE on the LHS. With this reinterpretation, we can determine the 
trajectory even in the interference region since the QSHJE, which is a partial differential 
eq., describes the interference; see ^13.31 



3. Application 

Based on the formalism developed in we explain how quantum phenomena relevant to 
a spinless massive particle appear. 



3.1. Energy quantization in 1-D space, Uncertainty relation, Tunneling, and Free 
particle 

Energy quantization in 1-D space The energy of the QSHJE with a confining potential 
is quantized to the same levels as that of QM, which is seen from ( 12.12b . 
Uncertainty relation between position and momentum We examine a particle in a 1- 
D potential well. Let the energy £ of a particle be £ > 0. Let the potential of the well 
be V{x) =0 for < L, V{x) » E for |x| > L. Let k = \/2mE/h. We assume that the 
particle has the lowest possible energy. Then we can take ^(x) ~ cos^x and ^^(x) ~ 
sin^jc as two linearly independent solutions of the corresponding SE. The momentum 
p = d§o I dx in the well is found using the formula (12.9b . Assuming I2 = 0, we obtain (note 



dx {^>0 + l2^>)2 + f^^>2 sin2g + /2cos2g l + (/2-l)cOs2g ' 

From which, assuming li > 0, Ap :— max. \p\ — min. \p\ is obtained as 

hn , 1 , 
^^^=A^I''-/;I' 

where the width 2L of the well was replaced with Ax. If /i ~ L2 or 0.85, the quantum 
mechanical uncertainty relation ApAx ~ /i is recovered. 

Tunneling From the energy formula ( 12.351) . we see that it can be q^ > even in a region 
where E — V{x) < 0, or that tunneling occurs respecting the energy conservation law and 
the particle picture. We note that we cannot determine the trajectory in the region where 
£ w y (x) since the semiclassicality assumption breaks down there. 

A free particle The trajectory of a free particle is not necessarily linear: q(f ) = vf (v = 
const.) since our E-L eq. ( 12.33b is of the 4th order. The deviation from q(f ) — vf, if any, 
is estimated to be very small, say, from Time-of-Flight experiments. When the deviation is 
small, the trajectory is approximated by q(r) = vf + Ct cos(a)/ + 6). Inserting this q(f ) into 
the energy formula ( 12.35b . we see that the energy E of the particle is expressed as 

E ~ imv^ + Ima^co^ . (3.1) 
2 2 
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Even if the particle is charged, the electromagnetic radiation does not occur Such ra- 
diation to occur, it has to be maP'O^ /2 > h(0. Substituting h(0 with (see footnote|3i, 
we see that a)^a^/2v^ > 1, which contradicts the semiclassicality assumption in which the 
oscillation is small. 



3.2. Energy quantization in 3-D space and angular momentum quantization 

Phenomena corresponding to the quantizations of I. (z-component of the angular momen- 
tum), (square of the angular momentum), and E (energy) of QM appear to the system of 
a particle in a 3-D attractive central force field V{r). 

The QSHJE of which is separated to three 1-D eqs. on the spherical polar coordinate 
system connected to the Cartesian one with x = r sin cos (p, y = rsinO sin (p, z = rcos 0. 
That is, the QSHJE 



f ^-^(VSo)^-VM+£ = 0, V(^ 
Zm K 2m m 



VSo)=0 (3.2) 



is, setting R — RipRgR, , So — So<p + §09 + Sor. and making use of the technique used to 
derive ( 12.61) from (12.5b . separated to 

4ot'- ^ 2m 4m^ ' 2m\ sin^ 2 4sin^0/ 

, 2io , /2m , , , , A \ 

-4;;^{^""'^>--2;;^(f^^-^(^»-;^)' ^'-'^ 

where ^rn^, are separation consts. in the first eq. of ( 13.21 ). The separation consts. in 
the second eq. of ( I3.2l i were set to zeros. 

First, we look for requirements for the solution of ( I3.2l i to be physically acceptable. 1) In 
1-D space, the solution So(^) of the QSHJE with a confining potential is required to satisfy 
( 12.12b . In 3-D space, if QSHJE is separable on the Cartesian coordinates to three 1-D eqs., 
we can express the 3-D version of ( 12.121 ). requiring the ( 12.121 ) on each direction, as 'The cor- 
responding 3-D SE has a L^-solution.' • • •(*) Since in this form, the requirement is indepen- 
dent from coordinate systems, we impose this requirement on the 3-D QSHJE with a con- 
fining potential regardless of the separability. 2) In 1-D space, the solution So(-*^) has to be a 
diffeomorphism from x to Sq. In 3-D space, this requirement is translated to: VSo 7^ 0, ±0° 
except at points at which the effective potential Vgff — h^(4m^ — cos^ 0)/(8msin^ 0) or 
Veff = y(r) +h^X/{2mr^) is singular (Veff — ±00). 3) In ^-direction, we require that i) 
d8iO(p/d(p is single-valued, ii) dSofp/ d(pd(p =nh/2 (n = ±1,±2,---), which is under- 
stood as follows. Consider the 3-D QSHJE having the confining potential separable on both 
Cartesian and spherical polar coordinates, that is the QSHJE with the harmonic potential. 
Then, on each of x-, y-, and z-directions, the ( 12.1 lb is imposed on the solution §0, (/ = x, y, z). 
Since §o{x,y,z) is a scaler, the integral value of the ( 12.111 ) on each direction is independent 
from the integral path taken in 3-D space. Taking the integral path of ( 12.1 lb for Sq, along 
the arc: r = const. >> 1, = 7r/2, n < <p < 2n, we see that the above requirement is 
equivalent to the ( 12.11b for §ox- If ^Soe /(30 7^ or d§,Q, /dr 7^ 0, « = or d§o^/d(p = is 
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allowed since VSq ^ for which. Summarizing; For a solution So of the 3-D QSHJE with 

an attractive central force field to be physically acceptable, 

★1) The 3-D SE corresponding to the 3-D QSHJE has a L^-solution. 

★2) VSo 7^ 0, ±oo except at singular points of effective potential. 

★3) d§0(p/d(p is single-valued. 

★4) jQd§o,p/d(pd(p = nh/2 (n = ±1,±2, • • •)■ If d§oe /dB or d$Qrldr 7^ 0, « = is 
allowed. 

Next, we examine for which values of m, A and E, the t^I) ^ 7^4) are satisfied. The 
solutions of ( 13.3b are constructed from the solutions of the 1-D SE's corresponding to (13.31) . 
which are obtained separating the 3-D SE corresponding to ( 13. 2t : 

_ VV+(V(r)-£)v/-0 (3.4) 
Im 

with \if{(p,9,r) = <l>((p)0(0)R(r). Explicitly, they are 



d(p 



2 



nr-^ = 0, (3.5a) 



1 d , . „t/0, 



• ^ (sin 0-777) + (^ ^)© = 0, (3.5b) 

--^ + (-^(£-y)--)R = 0. (3.5c) 



Allowed values for m If m 7^ 0, the general solution of ( I3.5ab is obtained from two linearly 
independent solutions <I> 1 = cos m 9 and <I>2 = sin m ^ as <I> = c 1 cos m ^ + C2 sin m (jO, where 
c\,C2 G C. The d§Ofp/d(p is constructed from which, using (12.15b . as 

d^Oq, _ h m(ciC2-CiC2) ^ ^ 

d(p i {ciCi+C2C2) + {ciCi—C2C2)cos2m(p + {ciC2 + C2Ci)sin2m(p 

Helped by the formula: f , . ,1" = ^=2 arctan ( i£d^}^B±^) , where > 

+ b^, we see that, to make Jq dSocp/dcpdcp — nh/2 true, it has to be m = ±1,±2, • • • . 
If m = 0, we take <t>i = 1 (see p|7]i and <I>2 = 9 as two hnearly independent solutions 
of ( I3.5ab . From which the general solution is obtained as <I> = ci 1 + C2<P- The d§o^/d(p is 
constructed using ( 12.15b as 

^§0y _ _^ (ciC2-ciC2)l ^ 

d(p 2i { ciCil^+{ciC2 + C2Cl)l(p + C2C2(p^) 

If ciC2 ^ c\C2, So(p does not satisfy the +3). If c\C2 = ciC2, Scp satisfies the *3) with 
d§o^/d(p = 0, which is allowed because 0-direction momentum can be non-zero as is 
seen below (see ^2) ). Thus, allowed values for m are m = 0, ± 1 , ±2, • • • . 
Allowed values for X and E Since the values of X and E for which the *1) is satisfied 
are found in standard textbooks of QM, we do not repeat the derivation of them (allowed 
values for X are X — + 1) (I = 0,1, - )^ those forE depend on V{r)). What we need 
to confirm is that the ^2) is satisfied for such A 's and £"s. We have seen d§Q^/d(p — for 
m = 0. To satisfy the *2), therefore, either of d§oe /dO or dSo,-/ dr has to be non-zero when 
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m = 0. We will see that, for any physically acceptable values of m. A, and E, (5§or/ dr ^0 
is not allowed, whereas (9Soe /dQ 7^ is allowed. 

For values of A and E which give a -solutions for ( 13.4b . the general solution of ( I3.5bl i 
and (13. Sel l are written as = C30con + C40div and R = csRcon +C6Rdiv, where 0con ■ ■ ■ Rdiv G 
K; the 0con (Rcon) is the 9- (r-) direction component of the L^-solution; the 0div (Rdiv) is the 
solution linearly independent from 0con (Rcon)- For the dSog/dO (d$Qjdr) to be non-zero, 
both of C3 and C4 (cj and cg) has to be non-zero (see (I2.15l l). However, Ci[{c(,) 7^ is not 
necessarily allowed because 3-D SE ( I3.4l i is not necessarily satisfied by the solution \// = 
<I>0R constructed from the 1-D solutions. Since r — Q and 9 —Q,Tt are singular points as is 
seen from the Jacobian J.J — r^ sin 0, the behavior of i/a = <I>0R there may be incompatible 
with the 3-D SE on the Cartesian coordinate system. We therefore examine the behavior of 
\lf — 4)0R there. 

As a preparation, we give the general solution of ( |3.5b| l around 0=0 for m = and 
A G C. It is obtained, using the Frobenius method (e.g., Ref.[T6]l. as0 

1 — cos 

0(cos0) =G0i(cos0)+//02(cos0) = GF(-/,Z + 1;1; ) 

, (-0,(/ + l), ^,^ 1-COS0 , '^'^ 1 1 2 1-COS0 

+^1^, (1)2 (i"^^+>^/^+7Tm-m^^(^^^ ' 

where G,H = const, e C; + = A; (•)() = 1, ==••••(• + r - 1) for r > 1; the 
Y^kloi' • • ) is understood as zero; the F{a,fi\y;z) is the hypergeometric function defined 
as F(a,/3;7;z) ir=o(«).(i3)rz7((7)r'-!) (kl < !)■ Note the l/{-l + k) = 1/0 and 
l/(/+l+A:) = 1/0 are cancelled away with a factor in the numerator of ^ttti-^. We 

know J + 1; 1; i^^f^) is a L^. function on < < ;/r for / e Z. Since Z = 0, 1, • • • and 
I — — \, —2, • • • give the same 0, we restrict the domain of Z to Z > 0. To list a few gen- 
eral solutions, 0(cos0) = G + //(ln i^^f^ - In iif^) for Z = 0, 0(cos0) =Gcos0 + 
//(cos In i^f^ - cos In i±f^ + 2 ( 1 - cos ) ) for Z = 1 . 

First, we show dSo^/dr = 0. § We assume that the potential takes the form: V{r) — 

fl_i/r + ao + fli'"H > where a_i,ao, • • • = const. G M. If Z = 0, ;i;(r) := rR(r) behaves like 

r or a const. 7^ for r 0, which is seen from the indicial eq. (e.g., Ref. [T6b. The xi^") 
which behaves like the latter is unacceptable. Assume the solution behaves like a const. — 
bo^Ofoi 0. Then, R behaves Uke bo/r for r ^ 0. The ^ = O0R = (1 + 0div)^o/'- 
(we ignore irrelevant coefficients; the explicit form of 0div is given above) constructed 
with which is not a solution of the 3-D SE ( I3.4l i as is seen as follows. If we insert y/ — 
(1 + 0div)^o / into (I3.4l i and integrate the eq. over a small sphere Q. of which center is the 
origin, we see the integral of the LHS equals 2h^bQ/m, or limn^o / {—h^'V^\j//2m +{V — 
E)\i/)dQ. ^ Ih^bo/m. Note jV^{l/r)dQ. = -An; jV^{&diybQ/r)dQ. = since 0div(f - 

'General solutions of )3.5bt for any m.A G C ai'e obtained from general solutions of the hypergeometric eq.: 
z(l - z)«" + (7 - (a + j3 + l))i(' - a/3i( = (e.g., Ref 17 1. General solutions of which for any a, j3,7 e C are 
found (exactly, were found on June 2010) in explicit forms with derivation on 
http://en.wikipedia.org/wiki/Frobenius_solution_to_the_hypergeometric_equation 

JThe following discussion is similar to that found in, e.g., chap. 9 of Ref. TJ] The difference is that ours involves 
©div and Rd IV- 
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9) - -0div(f + 0); limn^o/(V-£)v/t/n = lima^o/ ^^(1 +®div)'-^ sin 0^/6^/9 = 0. 
On the other hand, that of the RHS is zero. It therefore has to be bo — 0. Whereas the 
X's which behave like r for r — > are allowed. If Z 7^ 0, xi'') behaves like r'+' or r^' for 
r — > 0. Solutions which behave like r^' for r — > are not allowed, which is is seen in a 
similar manner as above. Solutions which behave like r'+' for r — ^ are allowed. Thus, for 
any acceptable values of /, only one from the two linearly independent solutions: Rcon('") 
and Rdiv('"), is acceptable. (The xi^")'^ which behave differently for r are linearly 
independent.) Accordingly, dSor/dr = 0. 

Second, we show that d§Qg/dO 7^ is allowed. To see Xj/ = 4>0Rcon = *J'(c30con + 
C4©div)Rcon with C4 7^ is allowed, we insert this y/ into ( 13.41 ) and integrate the eq. over 
small volume: < 6 < 5 = const. << 1 as in the case of the r-direction. Then, we see 
Xj/ = <t>(c3©con + C4©div)Rcon does not Contradict (13.4b . which is understood as a result of 
the slow divergence of ©div for — > 0, tt. Thus, we see (3Soe /dO 7^ is allowed. 
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3.3. Interference 

We determine, in the setting of FiglH the trajectories of electrons in the interference (i.f.) 
region, and the density distribution of them on the screen. We calculate the trajectory with 
(12.36b . the density distribution from the hitting positions. We will see that the density dis- 
tribution differs from that of QM by a few percent. 

• Apparatus We study with the apparatus of FiglT] Explicit values are: 

kx = ±4.99 X \Q^mm-\ ky = 0, k, 1.45 x \{fmm-\ 
filament diameter: 4.15 x lO^'* mm, 

aperture size: —0.2mm <x< 0.2mm, —0.2mm <y< 0.2mm, 
the distance from the filament to the observation screen: 33.77mm, 
where kx{ky,kz) is the wavenumber in the .a; (jjz) -direction after passing the electron 
biprism. The electron source is assumed to be located infinitely far from the biprism. With 
these values, we simulate the apparatus of Ref . [T] 

• Incoming beam In our formulation, only a superposed Gaussian beam (by 'beam', 
we mean the SB's solution) is allowed as a source side solution of the corresponding SE, 
which we call the incoming beam. The incoming beam is a solution of the free SE. The 
plane wave, the Bessel beam (e.g., Ref. TF') and the Gaussian beam (e.g., Ref. 19) are can- 
didates. The superposition of different energy solutions is not allowed since the QSHJE, the 
corresponding SE too, has a fixed energy; see S|4] From the plane wave, we cannot construct 
non-zero momentum in the .JC-direction behind the biprism (see (I3.18l l). The Bessel beam 
cannot be the incoming beam because the transverse amplitude of which changes sign as 
|r — ro|, where ro is the beam center and r is the transverse distance from ro, increases. We 
therefore take the Gaussian beam. The Gaussian beam, which is a solution of the paraxial 
form of the free SE, is written as 



-exp( j—\ )exp;fc(z- 



wiz) 



Vv2(z) 



[x-xpf + iy-yof 0(z). 
23?(z) k ■ 



(3.6) 



where wq is the beam width at the beam waist position, w is the beam width at z, (jcoj^o) 
is the beam center coordinates, k — {k^ + ky + k^), 3? is the radius of curvature of the 



source 





b i pr i sm 
aperture 



screen 



Fig. 1. Experimental setting: Electrons are emitted from the source at the left side end and accelerated toward 
the screen. Electrons passed through the upper (lower) slot change the moving direction downward (upward) due 
to the voltage applied to the central filament. 
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wavefront, and is the Gouy phase. The z-dependence of w, 3?, and is written using 
the Rayleigh range zr := kwl/2 as vp-{z) = w^(l + (z/z/j)^), 3?(z) = z(l + {iR/zf), and 
0(z) = tan^'(— z/zff); see, e.g., Ref. fT9]for details. We consider that the beam waist posi- 
tion, at which 3? = oo, is adjusted to z = with focusing lenses on the source side. Because 
preliminary calculation suggests that in the real world the transverse beam profile is close 
to rather than g^t'^^'o) , we introduce a superposed Gaussian beam. In the present 

calculation, we assume the simplest of which: two components Gaussian beam (see ( I3.15l l). 
We assume, although we don't exclude other possibilities (see pl24]i. that the incoming 
beam accompanies an electron at its center. The conclusion that the diffraction intensity 
distribution differs from that of QM is kept intact regardless of the transverse electron dis- 
tribution in the ensemble of identical incoming beams (see footnote[s]i. 

• Momentum field behind the biprism^ We determine the momentum field behind the 
biprism caused by the Gaussian beam (I3.6l l, from which the momentum field caused by the 
superposed Gaussian beam is easily obtained. 

The momentum field 
p(x) = VSo(x) behind 
the biprism is constructed 
from the solution 'Pdiff 
of the corresponding SE 
using (ITTSI l. The ^diff 
is the superposition of 
two beams, one which 
originates from the vir- 
tual source located at 
the lower left far from 
the filament and was 
diffracted by the knife- 
edge open downward at 
z = (no voltage applied), the other which is obvious from symmetry; see Fig|2]and Ref. 
l2l We use three coordinate systems (x,y,z), {x,y,z), and {x,y,z) which are related with 

IctX l^xZ ^ ^ kj^^x ~t~ k^z ~ k^x ~1~ kj^z ~ ~ kj^x ~\- k^z _ 

We consider two Gaussian beams 'I'g and '^q. The '^c (^g) propagates parallel to the z (z) 
axis toward the z+ (z+) direction. The '^c is diffracted by the knife-edge in the {x,y) plane 
at z = open downward from i = — d/2 to —°° (d is the diameter of the filament). The '^c 
is diffracted by the knife-edge in the {x,y) plane at z = open upward from x = +d/2 to 

-|-oo, 

''When the beam width wq is large enough to form a reasonable interference contrast on the screen, the electrons 
which passed close to the filament accumulate in a small area on the screen, which contradicts experiments. 
'Although we find the formula which gives the Gaussian beam field diffracted by the kiiife-edge in Ref, |20| it is 
too elaborate for our purpose. We therefore determine the diifracted field in a simpler setting. 



virtual source ^ 



virtual source 



^ / downward open 
^ knife-edge 




\ upward open 
\ kn i f e-edge 



Fig, 2, The coordinate systems 
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The disturbance U{P) of "Pg at the point P — (X,?,Z) behind the downward open 
knife-edge, where {X,Y ,Z) are coordinate values on the (i,3',z)-system, is given as (e.g., 
Ref.EDchap.S) 

m^-^ r'\l r.^exp(- (^"-^»)^ + (^-^»)V /..) , (3.8) 

As' J -co J-oo Wq 

where A: ampHtude of the incoming beam; A(= 2n/k): wave length; coordinates 



on the {x,y) plane at z = 0; = y (1 - |)2 + (? - fj)2 +Z2 , 5 ' = y/x^ + Y^+Z^ . We 
consider only the case: yo —0 and ? = 0. In this case, U{P) of ( 13. 8t is written, with the 



Fresnel approximation, as (see [Appendix 1 1 

U{P) 



/ k^l^ k^T^ k^l^ k^l^ 

X ((cos( — ) . M, +sin(_) . NO + /(cos(— ) • N, - sin(_) • M,) 

where l^X/S\ a = k{\ -F)/ {ns '), Mj = Mi(l,Z) and Ni = Ni(l,Z) are slowly 
changing (compared to cos,2kxx) functions of X and Z, explicit forms of which are given 
in [ Appendix II Rewriting U{P) as U{P) = ^^iE^^£l^i^(Mi + /N,)exp(-/|0-), we see 



that the field intensity \U (P) P caused by the lower-slot-passed beam is ^ (M^ + N [ ) (note 

In a same manner, U [P) is written as 

U{P) = e'^-' 
V2(l-P) 

/ A;2f2 k^l^ k^l^ k^l^ 

X ( (cos(— ^) • M2 + sin(— ^) • N2) + /(cos(— ^) • N2 - sin(— ^) • M2) 
V zTza Ina 2na 2na 

where what M2, N2, /, 0; represent are obvious. 

The diffracted field ^diff at P = (X,y,Z) behind the biprism is the superposition of 
t7(l,?,Z) and^(i,?,2): 

^diff = f/(i') + &(i')=e'*^''((Mi+/Ni)e-'^ + (M2 + /N2)e-'^) , (3.9) 

where j3 = {k^^)l{2na) and ^ = {kH ^) / {2Kh); the values of 1 - Z,l - 2 are given 
fromX -Z using (lOi.FI 

The momentum field VSq is obtained from ( 13. 9t as (there are implicit fixings of param- 
eters in the following formulas; see pl24li 

„5 Ti'Pd.ffV^diff-I'diffV^diff 

p = V§o = T ^ , (3.10) 

I 2>Pdiff1'diff 



m^g-m/4/^2(l-/2) ~ Ae-'''''^/V2 = const, is omitted. 

Note s ' = v^F+F+l2 = \/F+F+l2 = I ' = y = Vx2 + f2 + z2 . 
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We give explicit form of px ipy ~ 0,/?^ ~ hk^ is obvious). In 



aso n ^diff^x^diff - ^diff^.v^diff , , , 

Px = ^^ — = ■ (3.11) 

dx I 21'diff'Pdiff 

We have 

I'diff'Pdiif = + + M| + + 2(M 1 M2 + N 1 N2) cos 2;t;,x + 2 ( M 1 Nz - N 1 M2) sin Ik^x, 

(3.12) 

in which j3 — j3 = —Ik^x was used, and 

*d,ff5r^d,ff-^d,ff<3x*diff , , 

= lm(^diff^.t1'diff 

2i 

Ni(9^M2 + Mi(9iN2 — N2(3xMi + M2<3j:Ni ) cos 
+ (^-MidM2-^idx^2 + ^2dMi + N25xNi) sin2A:;,x , (3.13) 

in which -^{ks' — ^) ~ k^, -^{ks' — j3) ~ —k^ was used. Thus we see that 
p., ^/i-RHS of (I3j3|)/RHS of (I^TI)) . The forms of Mi,Ni,M2, and N2 at P = (X,0,Z) 
on the coordinate system {x,y,z) are found in [Appendix I| 

• Velocity field The VSq is connected to (q,q,q) with ( |2.36l l which is reduced to 
V§o = mq since q is large. Thus, the velocity field, which gives the particle's velocity at x 
if the particle comes to x, is given as 

1 d§o fi RHSofdHll . ^ , hk, 
m dx m RHS of p.l2p m 

• Calculation of the electron distribution on the screen We will calculate the electron 
distribution on the screen chasing trajectories obtained solving ( |3.14| i. 

We assume that the incoming field is the two components Gaussian beam (see the 
sentences below (I3.6l l): 

% - (0.8exp(- (--°);3;(r^°)^ +0.2exp(- (---»); + ;^-^»)^ ) 

X exp(fc) . (3.15) 

where the unit of length in (jc,}') -plane is pim. 

Since the Fresnel approximation breaks down in the vicinity of the biprism, we cannot 
follow the trajectory beyond the biprism. Therefore, 1) we have to assume some initial 
electron density distribution somewhere (our choice is z = Amm) in the Fresnel region, and 
2) we have to determine the diffracted field with which the electron located at {xi„i,yini) 
at z = 4mm interacts, or equivalently, the center of the incoming beam from which the 
diffracted field originates. Regarding the 1), we assume that the electron density at z = 4mm 
is equal to the field intensity formed by the uniform incoming Gaussian beam flux, which 
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turns out to be equal to that formed by the incoming plane wave up to a const, factor. 
Regarding the 2), we assume that the electron of which position is (jc/n,-,}',-,,,) at z — 4mm 
interacts with the field which originates the Gaussian beam came along {xq — Xini — j- ^ 
4mm, yo =yini). 

The electron density /lower (j^) on the screen caused by lower-slot-passed electrons is 
reconstructed from the trajectories as follows. Let the initial jc-coordinates of n sufficiently 
many electrons be JCi /,,,- < • • • < x,.„„ < • • • < x„.„„. The corresponding hitting x-coordinates 
be xi.i,i, < ■ ■ ■ < Xf.fjif < • • • < x„.f,ir- Let the field intensity at z = 4mm be |^diff(-x)P- Then 
/lower at x;./„Y is given as 

, / N l^diff {Xi.im)\^ {x{i+l)-ini-X{i-l).ini) , . , ,n , 

/lower ) °^ 7. 7^ — ^ — 7 (2 < i < n - 1 ) . (3.16) 

^Ixihit [X{i+l).hit—X{i-l)-hit) 

The /upper(^) by upper-slot-passed electrons is given as /upper(-^^) — /iower(^-'^)- 

Figgis the assumed electron density at z = 4mm. Fig|5]shows the calculated trajecto- 
ries of electronsH The upper graph of Fig|4]is the electron density pattern calculated from 
the trajectories. Ifl The middle graph is the field intensity calculated with the uniform incom- 
ing flux of two components Gaussian beam ( l3.15b FI The lower graph is the field intensity 
calculated from the plane wave incoming beam. Note that the envelope shapes are almost 
same in the middle and lower graphs. Note that the largest Fresnel peak of the electron 
density is shifted to the outerward at the arrowed positions. This is explained as follows. 
The electron runs almost linearly parallel to the geometrical border, whereas the Fresnel 
peak of the field runs alongi = const. y/i; see, e.g., Ref.|2T] §8.7.3, see Fig.|5j 

• Velocity in the x-direction behind the biprism As is seen from ( I3.14l i. the velocity 
in the x-direction in the i.f. region is slower in the both-slots-open case than that in the 
one-slot-open case (kx{M\ + Nj) — kx{M2 + N2) in the numerator is essential) because the 
oscillation energy is supplied from the x-component of the translational energy (see (13.1b ); 
note 3-D QSHJE separates into three 1-D eqs. behind the biprism where the electric poten- 
tial of the biprism disappears. 

" The density of the lower-slot-passed electrons at z = 4mm is M j + Nj from )3.12t . The M 1 and N 1 are given in 
jot . Note 

1°° ( j e-'-'^---'«iy(t)d'[fdxo = 1°° dxo j dyj rfz e-(-"»'7(y)e-("'»'V(z) 

= / f{y)dy j f{z)dz j°° e-'-''--'")' ivo =co«i/.( j fi'c)dtf . 

°It took about 20 min. to calculate one trajectory on a 3GHz class personal computer with Mathematica. 
P The upper graph was reconstructed from the trajectories of the lower-slot-passed electrons of which initial x- 
coordinates .v,„, 's range fromx,,,, = —2.00011m to — 0. 105/Jm, with 0.01/im step for a,,,, = — 2.0/im ~ —0.930^lm 
and with 0.005/im step for a',„, = — 0.925 ;Um ~ — 0.105 ;U»j. The upper slot data were obtained reversing the lower 
slot data. 

1 The field intensity I(x;xo) caused by the two components Gaussian beam of which center {x,y) is (ao,0) is ob- 
tained from j3.12) replacing Mi ~ N2 with 0.8M i +0.2Mii ~ O.8N2 +O.2N21, where what M n ■ ■ ■ N21 represent 
are obvious. The field intensity I{x) is obtained as /(.v) = J'i^I{.x;xo)d.xo, where c should be set sufficiently large. 
In the present case, we have set c = 6.4^im. 

' The lower graph was calculated using 13.12) setting wq appearing in Mi ~ N2 to t». 
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We estimate the mean velocity in the x-direction at the screen position from the fringe 
visibihty (f.v.) of the interference pattern. The f.v. is defined as f.v. — {Imax — Imw)/{Imax + 
Imin), in which I,mx {Imin) is the maximum (minimum) electron density on the screen. We 
consider that / |l/i| ~ m|(5So/(9x|^'. For simplicity we assume that, in the i.f. region, 
the ^diff is expressed as 

"Pdiff = e'^-'(7upe'*^'-^" + Jine-""^') , (3.17) 
in which 7up,-/dn = const. > 0. Then, dSo/dx is given as 

d§o_ hK{Jl^-Jl) ,,,,, 

dx J^^+J^^ + 2JupJdBCos2kxX 

From which we have f.v. = 2JnpJdn/{Jup+Jdn)- We define from dTTSI ) the mean velocity 
in the x-direction x,„„ at the screen position as mxmn — hk^ (7up ~ -^dn ) / (-^up + -^dn ) ■ Thus, we 
have 



f.v.^ )2, (3.19) 

V 

where xq — htx/m. Some values are, fori,„„/io = 0.999, f.v. = 0.045; for i,„„/io = 0.990, 
f.v. = 0.14; for Xmn/xo = 0.900, f.v. — 0.44. The f.v. calculated from the upper graph of 
Fig|4]is ss 25%. Accordingly, the mean velocity in the x-direction at the screen position 
is, from ( 13.19b . ~ 97% of that in the one-slot-open case. Thus, in the present case, the 
narrowing (compared to the one-slot-open case) of the envelope width of the i.f. pattern is 
not seen clearly. 

• Implicit assumption onA^D We calculated the trajectory assuming that an incom- 
ing beam accompanies a particle at its center We however have no reason, at present, to 
exclude the possibility that the particle is located off-centered. In such a case, it can hap- 
pen that the velocity field on the trajectory given with ( I3.14l i changes sign at some point in 
the i.f. region. In other words, it can happen that the particle changes the moving direction 
from upward (downward) to downward (upward), which is what happens in B.I. (see Ref. 
|2]§5.1). In our formulation this moving direction change is not allowed because the J-O 
momentum p^jj is conserved in the region where no classical potential exists; see ( IC.6b . 

We show how we avoid the moving direction change. Remind, in 1-D space, the mo- 
mentum pv = d§o/dx was determined, giving two linearly independent solutions ^ e M 
and e M of the corresponding SE and Mobius parameters A ^ D, as (see (12.8b ) 

^§o_/ h{AD-BC)W 

dx 2{A^>o+B^>){c^"^+D^>) ' ' ' 

Also, we saw that the 5So / dx determined from [a + /c)^^ + (/? + /d)^ = Re^^^l^ is equal 
to that determined from (13.20b with A = a + ic,B — b + id,C — a — ic,D — b — id (see p|7|i- 
We assume again for simplicity that the I'diff in the i.f. region is given with ( 13.171 ). Then, 
writing x-component of ^diff as (/up + ^down ) cos k^x + i (7up — ^down ) sin k^x, and identifying 
(7up + -/down) cos A:,vX and (/up — -/down) sinfe^x as ^ and of (13.20b respectively, we see that 
calculating the x-component p^ of the momentum using (13.18b is equivalent to choosing 
{A,B,C,D} — {/,!,— /,!} in ( 13.20b . If we choose {A,---} — {— /,1,/,1}, the sign of px 
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is reversed {px —fx)- We can thus avoid the velocity reversal for the lower-slot-passed 
electron choosing {A, • • • } = {/, 1,— 1} if 7up > -/down, or {A, • • • } = {— /, 1,/, 1} if < 
■/down- The choices for the upper-slot-passed electron are obvious. 
Accordingly, the formula ( I3.14l i is modified to 

^= — 5- = =^- pTTo fi-, ^ 3' = 0, z= — , (3.21) 
m ax m RHS of (|3.12l) m 

where the 'h- (-)' sign is for the particle which passed the lower (upper) slot. Note that the 
intervals on 1-D space (jc-direction) divided by the point at which px = are disconnected; 
the Px is never zero in our formulation as is seen from ( 12. 9t . The values of the global consts. 
A ^ D therefore need not be the same in each interval. 

• Experimental verification We saw above that our formulation makes different pre- 
dictions from QM.Hour formulation therefore testable. Examples of such predictions are 

i) When the f.v. is small, the x-coordinate xp on the screen of the largest Fresnel peak of the 
lower- (upper-) slot-passed electrons will increase (decrease) almost linearly as the distance 
Zsc between the screen and the biprism increases. While QM predicts that xp increases 
(decreases) non-linearly as Zsc increases; see the paragraph below (13.16b . see Fig|5] 

ii) The larger the f.v. is, the narrower the overall i.f. pattern width is, which follows from 
(|3.19t . While QM predicts the overall i.f. pattern width is irrelevant to the value of the f.v. 

We have encouraging data. On Fig|6] although the narrowing is not seen according to 
the small f.v. (w 0.02), the peak shifts seem to appear at the upper left and upper right 
corners. Similar peak shifts are seen also inRef.landRef.H] 



Even if the electron density in the (ensemble of the identical) incoming beam 1', is assumed to be proportional 
to 1 1*,, I ~ , the predictions below are still valid because the velocity in the ^.-direction is given not with J3.14t but 
with j3.2H . If the electron density at z = 4mm is given with |1'diff|~ and the velocity is given with B.14K the i.f. 
pattern (electron density) on the screen will be the same as the field intensity. Note if mq = VSq, from the second 
eq. of i2.4\ . = l^^diffP represents electron density. 
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4. Difference between Bohm's Interpretation and the Present Formulation 

We first examine how our formulation differs from B.I. , then show that the argument against 
introducing the trajectory to a relativistic spinless massive particle in B.I. does not apply to 
our formulation (introducing the trajectory to a relativistic Dirac particle in B.I. is possible; 
see, e.g., Ref.|2]l. 

B.I. extracts the particle from the solution of the time-independent SE of QM writing 
i/A = Re'^^l^ and interpreting VSo/m as the velocity of the particle: VSq = 'WQ- While our 
formulation, as a modified CM, has the concept of the particle from the outset. In our 
formulation, the trajectory is determined with the E-L eq. ( 12.331 ) or with the eq. ( 12.361 ). 

In B.I., the physical meaning of in = Re'^^l^^ is the same as that given in QM 
since B.I. is an interpretation of QM. Accordingly, represents the particle density (in 
an ensemble consisting of systems of the same y/). To make the particle number be time- 
independent in a time-independent potential, therefore, it has to be VSq = 'wq according 
to the second eq. of ( I2.4l l. While in our formulation the physical meaning of — Re'^°^'' 
is carried by the R and the Sq on the RHS. The §o is a kind of 'field' defined all over x 
regardless of the particle position. If a particle comes to x, V§o|x and (q,q,q)|x satisfy 
( |2.36t there. The R also is regarded as a kind of 'field' in our formulation (see i)3.3l ). At 
least, there is no reason for the R^ to be equal to the particle density. 

In B.I., the solution set {R, So) of i2.4\ is acceptable if in = 7?e'^''/s quantum me- 
chanically acceptable. The superposition of the different energy y/'s is therefore acceptable 
since it is acceptable in QM. While in our formulation, it is unacceptable. In our formula- 
tion, in accordance with CM, to a single action 'field' corresponds a single particle. In the 
case of a particle moving in a time-independent potential, the energy of the system is const. 
The superposition of different energy y/'s is therefore unphysical. Accordingly, the wave 
packet does not appear in our formulation, nor is it necessary. 

In B.I., once the \j/ is given, the VSq is uniquely determined from \j/ = Re'^'^l^\ In 
our formulation, even if is given, VSq is not uniquely determined. In 1-D space, VSq 
is determined giving and the Mobius parameters A ^ D; see ( 12.8b . For example, given 
gives uniquely VSq = hk. While in our formulation, the VSq can be even —hk 
depending on the values of A ^ D; see ( I3.21l i. 

In B.I., the velocity of the particle in a 1-D confining potential is zero (we restrict 
our discussion to 1-D space for simplicity). In B.I., the solution of the 1-D SE of the 
system has to be an L^-function since it has to be so in QM. Let e M be two linearly 

independent solutions of the SE. Let ^ be the L^-function. Then, is not an L^-function 
(see i on p|6]l. Thus, in the form: = I' = Re'^°^^\ the So is const, since ^' is real. From 
which q ~ VSo/m = follows (e.g., Ref.|2]§6.5.1). While in our formulation, it has to be 
\|/^a^^ + b^'^ with 0^a,beC (see ^TTT i. Accordingly, in the form: i// = /^e'^o/fi^ the So 
is not const. ElFrom this So, the velocity is obtained with the eq. ( 12.361 ). 

Lastly, we show that the argument against the relativistic trajectory does not apply to 
our formulation. In the spirit of B.I., the relativistic trajectory of a spinless massive parti- 
cle would be introduced as follows, i) Rewrite the unknown solution of the Klein-Gordon 
eq. (K-G eq.) as 4>(f ,x) = /?(f ,x)e''^(' '')/^'. ii) Insert it into the K-G eq. iii) The result- 
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ing eq. is ( |2.3l l. iv) The four-velocity is given as Mu^ = —[d^S-\- {e/c)A^^), where 
is given in ( I2.3al i. The f in ( I2.3bb is understood as the particle den- 
sity in an ensemble. The reasons why this method does not work are i) f is not positive 
definite, and ii) the velocity vector can be space-like; A space-like velocity is constructed 
using a superposition of the different energy solutions (e.g., Ref. 2 §12.1). In our formu- 
lation, however, i) f does not represent the particle density, ii) in the time-independent 
electromagnetic field, the superposition of different energy solutions is unphysical (see the 
second last paragraph). Thus, it seems that we need not be so pessimistic about the possi- 
bility of the relativistic trajectory of the particle, although the theory construction is left to 
the future. 

5. Concluding Remarks 

We have formulated a quantum theory of a spinless massive particle as a generalized (or 
modified) analytical mechanics, i) We modified the classical H-J eq. to what we call the 
QSHJE using the extended action of the diffeomorphism group, ii) We determined the 
semiclassical Lagrangian consistent with the QSHJE; we found that the Lagrangian is of the 
second-order: L(q,q,q) and that, within our theoretical framework, there is no Lagrangian 
beyond the semiclassical level, iii) Based on the QSHJE and the Lagrangian, we explained 
the appearances of quantum phenomena relevant to the particle keeping the particle picture 
throughout. 

Thus, one of the conceptual difficulties of quantum mechanics: the impossibility to 
maintain the particle picture, was resolved. 

Our formulation is testable (with a present-day technology) because the particle density 
in an ensemble differs from that of quantum mechanics. 
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Fig. 3. Assumed electron density atz = ^mm. 
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Fig. 4. Interference patterns on tlie screen; Upper graph: electron density reconstructed from the trajectories; 
Middle graph: field intensity caused by the two components Gaussian beam )3.15t ; Lower graph: field intensity 
caused by the plane wave beam; Note that the electron density is larger than that expected in QM (that is, larger 
than the field intensity) at the arrowed positions. The marks on vertical axes are arbitrary. 
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Fig. 6. Experimentally obtained electron density pattern and the calculated field intensity pattern reproduced 
from Ref. 7. 
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Appendix A. Extended action of Diff(RP') on J2(RIP') 

We give the mathematical background of the extended diffeomorphism action appeared 
in f|2T| In the following, we naively apply the results of the theory of finite-dimensional 
group extensions to the infinite-dimensional Lie group Diff(MP'). For details on finite- 
dimensional group extensions, see Ref. [m[T2l[T3l etc. 

Let the diffeomorphism group of 1-D projective space MP' be Diff(RP' ), the covariant 
tensor fields of order two (equivalently, the group of tensor densities of degree two) Q on 
KP> be J2(KP')- 

We consider the classical action T of Diff(MP' ) on J2(]RP' ) as a special case of more 
general action, in which T is given as T : Diff(MP^) x 3^2 (M^) J2(KP'), (y'\<^{y)) ^ 
Ty-i(p{y) = (p{y{x)){dy/dx)^. Let E be an extension of Diff(MP') by 5'2(KP') or ^ 
J2(KIP') ~> E ^ Diff(MP') ^ 1 be a short exact sequence. Using the set-theoretic bi- 
jection between E and 3^2(R1P') x Diff(RP^), we write an element of E as {^,g), where 
e J2(KP') and g e Diff(RP^). Let s{g) {dg,g) G E and := {(j) ,1) £ E, where s{g) 
is a section (a representative element of g in E) and is an injection from 9^2 (KP^) to 
E. For a given action T, the operation rule in E is given[j as {(j)i,g){^,h) = (0i + Tg(f)2 + 
f{g,h),gh), where f{g,h) is a factor set which satisfies s{g)s{h) = i{f{g,h))s{gh). As a 
diffeomorphism, we require the s{g) to satisfy s{g)s{h) = s{gh) or the above sequence to 
spHt; accordingly, f{g,h) = 0. In E, the action of Diff(RP') on J2(RF') is given as the 
product of s{g) and/(0): s{g)i{^) = {dg,g){^,l) = {dg + Tg(l),g), or T(ag,g)^ ^ dg + T^^. 
If we set dg = 0, the classical action is recovered. In the present case that (p G 9^2(KP'), 
we can set dg^^ — a{g{x),x}, where a — const. G R and {g{x),x} is the Schwarzian 
derivative, because it satisfies — dg^^ + T^-idh^^ required fromi(g^')i(/i^^) = 

s{g-^h-^). In this case, the action T^d^-i^^-i) of Diff(RP') on ^2(^1"') is given as ^21 

Tidg-K g-^)^iy) = idg/dxf^igix)) + a{gix),x}. 



'A covariant tensor field of order n on 1-D space M' is written as (j> = (j>{x){dx)", where x is a local coordinate. 
The field is also called a tensor density of degree n. Equivalently, a tensor density of degree n is defined as a 
section of the line bundle (T'JVt')®", where r*M' is the cotangent bundle (e.g., Ref. 10). 

"The group operation in 5'„(M') is written additively (e.g., 03 = 0i +02), that in Diff(]Vt') and E is written 
multiplicatively (e.g., ^3 = gigi)- Accordingly, the unit elements are written as 0, 1, and 1 respectively. 
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Appendix B. Derivation of the Relativistic Quantum H-J Equation 

We derive ( 12.3b . Our Minkowski metric is diag(l, — 1, — 1, — 1). We use the Einstein sum- 
mation convention. We work with the Gaussian system of electromagnetic units. 

Fkst, we examine the effect of the extended diffeomorphism action on the simpUfied 
relativistic H-J eq. We set in the relativistic H-J eq. of a charged particle:'^ 

(^mS + %)('5^S + = ^(^ - - 1a)2 = (B.l) 

c c at dx c 

to zeros, then restrict the eq. to 2-D time-space (x'^{= cf),x') to obtain 

We examine the effect of the extended action of a diffeomorphism y : i-> y{x^) upon 
the second term of the LHS of ( IB. 21 ) with the first term left intact. Note, with the classical 
(not-extended) action, ( IB. 21 ) is transformed to the eq. on the coordinate system {ct,y) as 

/l ^§(f,x'(y)) x2 dy 2( d§it,x\y)) Y 
\c dt ) ^dx^' \ dy ) 

How is (IB. 2b transformed if the action is given with ( 12.1b ? Let cf' — > (j''' — > <7'(=x') be a 
transformation among coordinate systems. By abuse of language, we use the same symbols 
^"j^'', • • • as labels of coordinate systems, as transformation functions, and as coordinate 
values. The coordinate transformation on 1-D space from q'^ to q'^ is given by the action 
of {dq' -\q' e Ext. Diff(Ml) upon {dS/dq' f G JiiM^) as (see dO) ) 

where a — const, e M. Successive transformation from q'' to q" takes ( IB .4b to 



m^c^ . (B.3) 



T, 



Using the cocycle condition: {q'^,q"} — (dq'' /dq")^{q'^ ,q''} + {q'',q"}, we see that (IB.5I) 
is equal to 

Ti,,^ -K," -^)Tid,^ -Kef -.)(^)' = (^)' + «{?^^"} • 
Thus, for the coordinate transformation q" — > q'* — > the eq. ( IB.2b on q'^: 

is transformed to the eq. on q'' as 
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Then, (IB.7I I is transformed to the eq. on q" as 

While the eq. on q"^' is considered as the transformed eq. from that on q*. Thus, the eqs. on 
q'', q'\ and q" take the same form: 



However, an equation upon the coordinate system q'^{q'',q") cannot contain a transfor- 
mation function q* to another coordinate q*. The q* therefore has to be replaced with 
some quantity defined on q'^{q'\q"). In ( IB.9al ). if q* =Aq'^ + B (for which {q*,q'^} = 0), 
§ = —cpot + piq^, where po = const. > and pi = const. 7^ 0, is a solution, which hints at 
replacing q* with §. This replacement is supported by the following statement: A physical 
quantity defined on 2-D Minkowski time-space {ct,q'^) which is relevant to a free particle 
and linear in ^' is almost singled out as S. Accordingly, we set q* = S.0Thus, if the action 
is given with (12.1b . (IB. 2b is transformed to 

When an electromagnetic field exists, ( IB.9ab becomes 



With the same reasoning as that in the last paragraph, diq* is singled out as di§ — ^A^ . We 
therefore replace diq* with(9iS — fA'. Thus, (IB.l lb is written as 



c dt c ' V'^jci c ' ^ (^iS-Mi) 2^ (^iS-Mi) ; 

(B.12) 

Next work is to determine the 4-D relativistic Quantum H-J eq. 

First, we will find the eq. valid for the electromagnetic potential A's which take the 
form: A = (A''(x''), A^(x'), A^(x^), A^(x^)). Let a set of all electi-omagnetic potential A's 
which take such forms on a Lorentz frame x be ^Isep- If A e A^^^, dBJj separates into four 
1-D eqs. with S = SoC^") + Si (x^) + §2(-^^) + S3(-f^) on the frame. □ On the assumption 



""Because q* e Diff(M'), as a result of equating S with g*, S (seen as a function of q"^ {q'' ,cf ) at any fixed time t) 
is required to be, at least, monotonic. 

"The idea to use A = (A" (a" ) , ■ ■ ■ ,A^(x^)) was borrowed from Ref . i23i in which Bertoldi et al. derived the 4-D 
relativistic quantum H-J eq. from 'the Equivalence Principle' elaborately explained in Ref.[3] 
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that the extended diffeomorphism group acts upon 9^2 each 1-D eq. is modified to the 
form similar to ( |B.12t . Collecting each modified eqs., we have 

where Oo, • • • , as G M and 

di{d,$ + ^A^) 3 ao(ao§ + fA») 2 _^ df{di2>-^-A') 3 a,(a,§-fA') 2 

^ (^oS + fAO) 2^ (<3oS + fAO) ^ ' ^ (^/S-fAO 2^ (a,S-fA') ^ 

(/ = 1,2,3; take no sum). Interestingly, say, (likewise for 2^,0^,2^) is written as = 
2did^Ri/Ri with the solution 7? i of (9i (7?^ (<9i S i - eA '/c) ) = 0. Then, considering Oq = 
CCi = fp'/l is realized in nature (the resulting eq. becomes Lorentz invariant if Oo = a,-, 
reproduces experimental data if Oq = a, = we sum up non-classical terms Q : — 

cxoQ° - UiQ in (IbTjT i to -h^d^d^R/R, where R = Ro{x^)Riix^ )R2{x^)R3{x^). The (I03b 
is therefore rewritten in the form: 

(a^S + -A^){d^$ + ^^'') = + ^^^^ (B.14a) 

^m/=0' / = _^(aM§ + £AM). (B.14b) 
m c 

The (IB.14I I is a system of equations of which unknown functions are S and R. Although 
(IB. 14b looks Lorentz invariant, the invariance is not assured for A ^ Tlsep since ( IB. 141 ) is 
nothing but a rewriting of ( IB. 13b . 

Next, we will find the eq. set valid not only for A e Tlsep but for arbitrary A's, that is, 
for A'' = A''(x'',x',x^,x^). We assume that such an eq. set takes the form: 

(a^S + -A,){d^§, + -A'') = m^c^ + tMUl +g (B.15a) 
c c R 

di^f^g. /^-^(^Ms + V), (B.l5b) 
m c 

where g{g) is a Lorentz scaler which vanish {g{g) = 0) at least when A e ^Isep- The g{g) de- 
pends on some of the Lorentz scalers constructed from (9^S, A, R and their derivatives. The 
eqs. to be gauge invariant, the A appears in the form: (^''S + -A^). Thus, g (likewise for g) 
looks like g{x) =g{df,{d>'§+^A^'){x),{d§ + fAf{x),R{xidf,d^R{x),d^Rd^R{x),---). 
Because g{g) = for A G ^sep (including A = 0), and the functional form of g{g) is com- 
mon to A G j4.sep and A ^ ^sep, ^(^) = for arbitrary A. 

Thus, we see (IB. 14b as the modified H-J eq. on the assumption that the extended dif- 
feomorphism group acts on each term of the LHS of ( IB. lb and call ( IB. 14b the relativistic 
Quantum H-J equation. 

"From d, {R]idiSi - fA')) = 0, if 5iSi - fA' ^0,Ri =±const. {d,§i - fA')-'/^ Therefore 
Ri ~ Ri ~ 2\. 5iSi-^Ai 2^ 3iSi-^Ai ')' 
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Appendix C. Higher-Order Analytical Mechanics 

We summarize the 3rd order analytical mechanics, in which the Lagrangian is a function of 
q, q, q, and q. Although it is possible to consider still higher-order analytical mechanics, 
we don't need them. Materials in this section are found in Ref l5]and|6l References l24land 
|25]are useful as supporting readings. See them for details. 

• Lagrangian fonnalisiTi^ As in CM, action S is defined as S := /^f L(q,q,q,q)t/f. Let 
5V be the variation of the curve q(f ). Then, we have 

5S= / L(q + 5V,q + 5V,q + 5 V,4 + 5'^)<ir- L{q,q,q,q)dt 



dq dq dq dq J dq dt dq dt-^ dq dt^ dq 



,dL d dL d^dL^^ ,dL ddL^^.,u ,<9L,,.. ,„ 



b 



From the first term in the RHS of ( IC.lb . we obtain the E-L equation: 

dL d dL d^ dL dL 



. — . (C.2) 
dq dt dq dt^ dq dt^ dq 



We define P(i), P(2), and P(3) as ( See (IC.ll l ) 

dL d dL d^ dL dL d dL dL 

3 

The P(i), P(2), and P(3) are called Jacobi-Ostrogradski (J-O) momenta. We define £c {L) as 

3 

£c (L) =P(i)q + P(2)q + P(3)q-i 

,dL d dL d^ dL, ,dL d dL, dL . _ 
-^^-diM^dt^^'^^^-didt^^^^-'^- ^^-"^ 

3 

The £c (L), which is reduced to {dL/ dq)q — L if L — L(q,q), is const, along the solution 
curve of the E-L eq. ( IC.2l i as is seen below: 

d 3 d , ■ , 

di " ^(P(i)q+P(2)q+P(3)q-'L) 

_,ddL d^ dL d^ dL dL d dL d^ dL 
~ dt dq dt^ dq dt^ dii dq dt dq dt^ d^ 

,ddL d^ dL, ,dL ddL,. ,ddL,. ,dL,.. 

dt dq dt^ dq dq dt dq dt dq dq 
dL dL .dL ..dL (dL d dL d^ dL d^ dL 



"^dq **aq **aq ^d^ Uq dt dq^ dt-^ dq dt^ d'^) °" ^^'^^ 

3 

We use also the symbol E for Zc (L). Along the solution curve of iC.2i , we have 

dp(i) dL d ,dL d dL d^ dL, dL dL d dL d^ dL d^ dL ^ 
— Li. = — ( 1 _) = \ \ _ — . 

dt dq dt dq dt dq dt^ dq dq dq dt dq dt^ dq dt^ dq 

(C.6) 
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If the Lagrangian does not contain q exphcitly, therefore, d'p^^-^/dt = 0. 

m (m) (0) 

• Hamiltonian formalism'^ We use a symbol q for d"'q/dt'" (m = 0,1,- ■ ■); q = 

q. If the Lagrangian is regular (|(92L/(3^'a^-'| ^ 0), Leg : T^M ^ T* (T^M), {q,''q ,■■ ■ /q 
) (Q7QiQ:P(i)iP(2)iP(3)) is a local diffeomorphism called the Legendre transformation; 
where T^M is the 5th order tangent bundle and T* [T^M) is the cotangent bundle on T^M. 
Note the 6th order E-L eq. obtained from the 3rd order regular Lagrangian determines a 

3 

vector field on T^M. Let H {L) o Leg . Comparing two expressions 

J"-i) dH ^ dL^ (n) ^ , dL.An) 

''"^^TT)'' *1 +a^'^PW' dH = -^dq+ qdp^„^ + {p^„^-—)d q 
d q d q 



1,2,3). 



(« 1,2,3), we have, using (|C2l ) and ( [C3] ). 

t/*"q'' _ dH t/P(„) _ dH 



dt ap(„) ' dt ^ (n^i) 



Appendix D. The q may depend on the integration parameters of the E-L eq. 

The formulas ( 12.201 ). ( 12.21b and ( 12.22b take the same forms in both cases: i) cj^'s (k = 
l,2a,2b,3,- ■ ■) are consts. ii) q's depend on the start point (.x:o(fo),io(fo)) and the end 
point (x(fi),x(fi)) of the action integral. Note the dependence on which is equivalent to 
the dependence on the initial conditions {q{to),--- ,q{to)) of the 1-D E-L eq. or on the 
{v,a\,bi,qo{:= initial position)} of ( 12.26b since any of them is a set of four integration 
parameters of the 4th order E-L eq. 

If q's depend on {xo{to),xo{to)) and (x(fi),i(fi)), the dClb becomes, for a 1-D La- 
grangian L(q,q,q;ck{x,x)) (we use a 2nd order Lagrangian for simplicity; we fix the start 
point (jco,io)), 

i<t\)Ati) . .. f4ti),Ati) 

d§ = L{q + dV,q+ dV,q + dV;cii{x,x))dt — / L{q,q,q;cii{x,x))dt 



)■*('! ) , d . d .. d 

{{8V— + 5V— + 5V—)L{- ■ ■ ■Mx,x))dt 
aq aq aq 

4n)An) d d d d^ d , , , 

aq dt aq dt^ aq 

+^^^±^)l[...-,c,{x+5x,x))5vymah) + ^J±::j^^ 

aq dt aq dq 

where 5x ~ 5y|(^^), 5x = 5y|(, <). Because the end point (x,x) is fixed in the E-L eq. 
determination process, the resulting E-L eq. is the same in both cases. That has the 
same form in both cases is seen from 

P^'^-^^O^^^l'^O^Tq-JtTq^''^---'''^'' + 
Similarly for p(2)- 
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Appendix E. Characteristics of the solution of the 1-D E-L Equation 

We examine the solution of the 1-D E-L eq. (12.20b . 

• The velocity never vanishes The qofs. solution of the E-L eq. ( 12.20b never take value 
zero unless ^ = as is shown as follows. Let q'(f ) be a solution of the E-L eq. Then, 

„ ^^4_V^^^V K^-^W-^ , (3fc-2)(fc-l)^^ ^ 
1 2 v' Q^V k{k-l)qq''-^ {3k-l){k-l)q'' 



where C2 = 3c2a + C2f„ are conserved (cf. (12.2 lab and ( 12.22b ). Removing q, we have 

Assume ^(f) for ^(f) 7^ 0. For <7 — > 0, no term of the RHS of (IE. 3b approaches ±00; 
otherwise, the LHS cannot be kept finite because the convergence speed of each term of the 
RHS is different. Thus, |^*^/<?^*^^| < const, for ^ — > 0. From which we see q approaches 
zero faster than or equal to q^. Consequently the q never take value zero. 

• The velocity is written as a Fourier series We show, for any semiclassical initial 
conditions (pfToll. the solution of the E-L eq. ( 12.201 ) of which C2 = 3c2„ + C2b < is writ- 
ten as q{t) = vf + 52~^|A„cos(«a)f + 0„). Since initial conditions are semiclassical, we 
have |^(fo)| << m\q{t(i)\^ / {h^ -Ici) and |^(fo)| << m^\q{ti))\^ / {-Ih^ci) by definition. 
Accordingly, all terms in the RHS of ( 12.20b except mq and —2c2h^q/{mcf) are small com- 
pared to the mq. The ( 12.20b is therefore written approximately as = mq — Ic-^qj {mq). 
From an approximate solution q(t) = vt +Acos((»f + 0) of it, we see that q takes value 
zero repeatedly and at time f at which q = 0, q (otherwise q{t) = vt is the solu- 
tion). Inserting ^ = into ( IE. lb and ( IE. 2b . and solving resulting eqs. for q and q, we have 
q = ± Jp^^-lmE )/mand q = ±{p(\) ± Jpf^^-2mE " 2m£'/ {2c2h^m^). 



Since and E are consts., at times at which q = + y p^jj — 2m£' ) /m, the values of 

q, q{— 0), and q are equal respectively, which means the oscillation is periodic and the a) 
is the fundamental angular frequency (freq.). Thus, the q{t) is written as a Fourier series. 

• The inf. order of An is (//v)" Since the a is the fundamental angular freq., 
Acos(djf + 0) =Aicos(cof + 0i) and (oAj/v^CoAi/v, - ■ ■ e 0(> f/v). The lowest inf. or- 
der oscillating term in (//v)2-terms: 2c2h^ -Sqq/^mq^) and ~2ci,h^{3q^ + 3qq)/{m^q'') of 
( 12.20b is const. ■ Ajcos(2m +29i + ^) generated from A; (B^cos(a)f + 0i)sin((Bf + 
Aia)^sin^(a)r + 0i) and Al(o^cos^{(ot + 9i). It sums up to zero with the 2a) term: 
const. •A2cos(2of + 62) coming from and —2c2h^q/mq'^. Thus, we see A2 G 0((//v)^). 
In a similar manner, we see A„ e 0((//v)"). 
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Appendix F. Determination of C2h in the 1-D Semiclassical Lagrangian 

We determine C2b in the 1-D Lagrangian ( |2.19t of which c\ = C2a — 0. We also determine 
V, fli, b\ and © of the trajectory ( 12.261 1 as functions of £\, £2, and k of Pqhj (12.28b . They 
are determined as follows. From the E-L eq., we obtain (O = mv^ j (h.\J—2c2h)- It gives 
C2h = —mv^/{4-E) because CO = kv (see next paragraph). The v in C2b — —mv^/{4-E) is 
determined as a function of ai, bi and E from (12.221 1. Lastly the ai and bi are determined 
from the comparison of pqhj and pel- 

• Determination of CO We determine the relation between k — \/2mE/h and CO com- 
paring PQHJ ( 12.28b and pel ( 12.23b . We expand pqhj of ( 12.28b as 



Pqhj ^ \/2m£(l -eicos2fex-e2sin2A:x+0(£^) + 0(e|) + 0(eie2)) 

= V2m£(l -eicos2/tx-e2sin2ytx+0(e^)) . (F.l) 

Because we are in the semiclassical world, £1,62 << 1; we wrote the O(e^) + 0(e2) + 

0(eie2) as 0(£^) with e = ^ el + , which is allowed since e^^, e|, |eie2| < (ef + e|). 
We expand pel of (12.23b . inserting c\ = C2a = and q = vt + f, as 

PEL =Pw + ^ + Oih') . (F.2) 

We note that O(ri^) of (|F2li is equal to 0{{f/v f). We know, among the oscillating terms in 
(12.26b . Ai cos{cot + 9i ) has the lowest inf. order with respect to (w.rt.) f/v ( [Appendix E 1 



The angular freq. of the oscillation of pel having the lowest inf. order w.r.t. f/vis therefore 
2co generated from A\cos^{cot + 9i). While the wave number of the oscillation of pqhj 
having the lowest inf. order w.r.t. e is 2k. We therefore have CO = kv. In passing, we see that 

0(e) = 0((//v)2). 

• Determination of C2b We insert ci — C2a —0 and q{t) of (12.26b into the E-L eq. 
(12.20b . Then, we have 

+ ( 4 ) + (I p) = (mo)^ + ^'^h^Y )A 1 cos (a)f + 01 ) + ( I ^ p )= . 

m V mv^ V 

From which we have 

2 4 2 

—m V —mv 

• Determination of p(^\) i^nd E as functions of a\, b\, CO and v Using // = CO^{{a\ + 
b\)/2- aibx5m2cot + {b\~ a\)coii2cot /2) + 0{5^), and so forth, where 8 //v, we 
extract non-oscillating terms from and E: 

. 2c2bh^ k 4c2fo/i^ if- , 3. , 

P(i)^m — + -T + 0{n )=m[v + f) 

^ ' m cp m 
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E = -mq^ ~ ^ + ^ + (^2^ ) = -m(y2 + 2v/ + //) 

2 m m cf 2 



_?«i!,/_W^6^^.^^,^5«^,/!_4^^ ^^^^^^ (H5, 
m m 

at the inf. orders up to 5^ to have (the h^- and /i'^-terms of (12.21 al l and (12.221) do not generate 
non-oscillating terms) 

2c2bn^ -4 {-CO^){aj + bl) ^ ^C2hh^ (0\a\ + b\) ^ , 
m 2 mv^ 2 

mv^ 2 ^ ' 



-mv +-m ( — ) h 4 hO(d ) 



m 9 mco (ai: + b,) h co (a,+b\i, , 
2 2 2 mv* 2 

• 77je ai,^7i as functions of e\, £2, ft) anii v To express ai,b\ of (IF.6I 1 and ( IF. 7b as 
functions of Ei, £2 and E (ai — ai{ei,e2,E),bi — bi{ei,e2,E)), we once again compare 
PQHJ and We compare (IF.2l i written in the form: 



P£L = 7^(1) + (ftjyfl? + /72cos(ft)f + +0(5^) 

= P{\) + —^^ — ' ' +ai/7isin2ft)f+ ' ^ cos2ft)f (F-8) 

^ ' mv^ 2 2 

with dF II ) written in the form 



PQfjj = v2mE[l — £icos2fex-£2sin2fex + 0(£ )) 
= V2mE{l - £1 cos2A:(vf +/) - £2 sin2fe(vf + /) + 0{e'^)) = \/2mE x 
(1 - £1 cos2ft)f - £2 sin2ft)f - 2kf{e2 cos2ftjf - £1 sin2ft)f ) + 0{{kf)^e) + Oie'^)) ■ (F9) 

to have 

2c2bh^C0'^ 



-aibi + 0{8^) = -V2^ {e2 + 0{e^)) (FlOa) 



2c2bfi'^0}^ a\-b 



-0{8'^) = -V2niE {ei + 0{e^)) . (FlOb) 



Note kf ~ //v << 1, which is seen as follows. The oscillation amplitude / is roughly 
given as increasing rate / times half-cycle time n/{2(o): / ~ f {2(o). Using ft) — kv, 
we see 2kf/ Tt^f/v. Note there is no term having angular freq. 2ftj at the inf. order 5^ in 
PQHJ of (|E3 nor in pel of dESJ. We therefore converted (5^ ) in (lF8]l to (S'*) in ( iFlOb . 
From ( IF. 10b we have (double sign corresponds) 



„5. 



V2mE , T p,/ 2\\ ,2 mv^V2m£ 



£l±^£2 + £2 +0(£2)),fc2^ ^ ^ (£i±^£2 + £| +0(£2)) 
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• Summary up to this point Thus, we have obtained (remind e = \ ei + ei) 



2 4 2 



(0 — kv, coh ~ — T — T = — ; (F- 11a) 



p = + 1^ + 0(5^) (F. 1 Ic) 

^ ' mv^ 2 

^ ^^^^^^ 

1 2 2 mv^ 2 

• Expressing C2b, v, CO, a\, b\, and p(i') as functions of£i,£2, and E We express each 
term in jF.l 11 1 as a function of £1,62, and E. We assume v > 0. Then the double sign in 
(IF. llbl) is fixed to minus. Thus, 

(£i+£ + 0(£')), bj= (-£!+£ + 0(£^)). (F.12) 

mv\/ 2mE mv\' 2mE 

Inserting ( IF. 121 ) into ( IF.l Idl ). and solving for v, that is solving mv^/2 + \/2mE{e + 0(e^))v 
— £■ = 0, we have 



v=^^(A/l + (e + 0(e2))2-(e + 0(e2))) = ^^(i - £ + 0(e2)) . (p.B) 

V ' m. 



m 

From ( IF.13I ). we have 

C2b^^^ = --{\-2£ + 0{£^)) , (0 = kv=—{\-£ + 0{£^)) (F14) 



V2;^(l+e + 0(e^)) . (F15) 
Inserting v of dF 13t into a\,b\ of (|F.12| i. we have 

^? = -^(1 +^)(^^ " ^) = "2^^^'^^ - £ + 0(e')) . 
However, it has to be (see (IF.l Oat ) 



—£2\/2mE = 7 — aib\ = (1 + £)a\b\ 

mv^ 2mE 



Therefore, if £2 > (double sign corresponds) 



ai=±]l- — {~£,~£) + 0{£^-), Z7i = ±W-^(ei-e) + 0(ei). (F16) 
If £2 < (double sign corresponds) 
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Appendix G. Inconsistency between Pqhj and p^i at (//v) 

To see the inconsistency, we compare the pqhj and pel at e^/^ — (//v)^. We have 

PQHJ = T—, ^^TT^-TT TV- = - e{\ -%)co^{2kx - 0o)) 

(sin;tji: + £2Cos;tx)2 + (l+ei)2cos2ytx ^ ^ 2^ ^ ^' ' 

+0(e2) = hk{l - e(l - ^)co&{2kvt + 2kf- Qq)) + 0{e^) 

= M(l - e(l - y )cos(2a)f - eg) - 2eyt/sin(2a)f - Qq)) + 0{e^) , 



where the 9q satisfies cos Oq = (e /2 + £\)l \ + ei)^ + & 



and sin Qq — e^I^ (£^/2 + ei)^ + £| ; remind e~^e1 + £|. While we have 

, 2c2hh\ p 5fp 3c,n' p / 4 

PEL^P(l) + 6~) + 2 8 - • 

The 3 co-term at e^/^ of pg//y is 



mEk^A^ 



-hkekAisin{3m + Bx - Bq) = ^^sin(3a)f + 30i) , (G.l) 

where we used e = mEA\/h^ (( IF.I6I 1 or ( IF.17l i) and 0g = —2Q\ obtained from the compar- 
ison of Pqhj and pel at — {f/v)^. While the 3a)-terms at (//v)^ of pel is 

— y!-^ cos(3cof + 01 + 62) + —r-r- sm{3m + 30i) 

- cos{3 cot + 39i) = Mix 



(4M2 cos(3«f + 01 + 62) + -rk^A] sm{3 cot + 30i ) + —^k^A\ cos(3«f + 30i 
V 4 4mv 



(G.2) 



where we used o = A:v, 2 0i = — Oq and C2b = — l/2 + 0(e). Inserting the E-L eq.'s solution 
( 12.27b with C2 in which replaced with C2b, we rewrite ( IG.2b as 

dal ^ k^A](i ^ + - — cos(3a)f + 30i) + — sin(3«f + 30i)) . (G.3) 

mv {-2c2hr- 4mv ^ 4 

We cannot make both of the amplitude and the phase of iG.3i to be equal to those of dG.lb 
adjusting C3. In a similar manner, we can show that OJ-terms at (//v)^ of pel cannot be 
made equal to o-terms at e^/^ of pqhj, though we don't work it out here. 
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Appendix H. Determination of of the 3-D free Lagrangian 

We determine of the 3-D Lagrangian L = + c^Tr'ic^ /mi^ + c^^Ti'^c^ /mr'c^^ + 

0{h^). We set = (pHHl. 

First, we prepare expressions of the E-L eq., p^^ and E obtained from the Lagrangian 
above (see jOl ). ( |236l ) and ( |C4] l): 



^2pf_( 2^1' , 16qq^; 8q4^; ^ 24q2(qq)q',- 
?M V q"^ q" q" q" 



'^1i + — \--A~^—6 Tg 



(4q2 + 8q4)^, 48(qq)^^,^ ^ ^ ^ ^^ ^^ 



qfi q** 



h^cf,r 2q4 3q2 8(qq)^ 
' q4 q4 ^ q6 

7q^ 8(qq)(q4) 4q^(q4) , 40q^(qq)' \ , 
+ ^3 V qio qio qio + qi2 ^"-^^ 

We determine Cj^. First, we see Cj^ < 0; otherwise a wavy Hnear trajectory cannot be a 
solution. For the negative and a semiclassical initial conditions, the solution of (IH.lb is 
written as g',(f ) = v,f +An cos(tO? + 9ii)+gi{t) • • • where v,- = const. > 0, ^; G 0(5^), 
5 := col Ai |/|v| << 1, I Ai I = ^/(Ajj +A21 +Aji), likewise for |v|. From the truncated E- 
L eq.: mqi - 2h^clf^i/{mq^) = • • • (*2), we have CO = mv2/(riv/(-2c^f )). Note gi{t) ^ 
£jJ""2A,„cos(n(Bf + 0,„) in general. Indeed, if we insert = y,f + A,„cos(na)f + 0,„) 
into (IH.ll i. non-oscillating terms appearing at the inf. order 5^ does not vanish in general. 
Inserting co = k\ ~ Y]=i Viy/TmEi/h into the (O above, we have 

m^v"* _ -m^{v\ + vl + v^f 



„3D _ _ _ 






To express v, as the function of £1,, £2,, and Et of (12.311 1. first, we write A,i as a function 
of (£1, , £2,). Then, we write v, as a function of A,i and Ei. We insert (* 1 ) into ( IH.2l i to have 

qi ni ifqi 

= const. + ^ ( -^(— -+anbn sinlcot + ^ cos2a)f + 5^) , 

m v^v; V 2 2 / 

where an — An cos On and Z?,! = —An sin Compaiing it with pgnji (cf. ( |F.9t ): 



Pg//// = \j2mEi{\ -£i,cos2a)f-£2,sin2a)f + 0(e2)) 
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at the inf. order 5 = £ ' , where £ = ^^^j ^J {efj + el,), we have 



from which we have 



mOJ^A^j 3 mCO^{aji+bfi) 



, , 0(5') = ^ 'L-^^' + 0(5') = JeI + e|. . (H.5) 

Next, inserting (*1) into (IH.3b . we extract the non-oscillating terms at the inf. order 5^. 
The result is 

E = imv2 + \m(O^K\ - ^^^f + 0(5^) = \my^ + \m(O^A\ + 0(5^) . (H.6) 



From ( IH.5l l and (IH.6b . we have 



v, = ^(v/r74^-y^) + 0(5^) = ^(l-y^) + 0(5^). 

(H.7) 



Inserting (IH.7b into ( lH.4b . we have 



1/ 2lf^,£,Je2. + e2 



^3^^__ - I" )+0(£t) 
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Appendix I. Explicit Forms of M i , N i , M2, and N2 

We rewrite (13.8b to an easy-to-handle form. We expand s in (13.8b with a power series of ^ 
and 77, then, takes the terms up to the second order (Fresnel approximation) to have 



s~s' -{ll +mf\) + 



(1 - I 2)^2 _^ (1 _ ,^^2) _ 2lmB,f] 



where I — X/S' , m — Y/s'. We consider only the case m = Y /s ' = Q. In this case, ( 13.81) 
becomes 



iA r^/^ ~ 



(^-ior + (77-350) 



,-1 ^2 



iks) 



iA 



-d/2 , /-+. 



(1-P)|2 



.1 +'-^H^ + ^^)) 



''0 

X exp(- 



We transform (IDI i to a tractable form. Setting | = ^/y/Si + kl/ (na), a^k{l- P-)/(ns'), 
we rewrite the real and imaginary parts of the ^ -direction integral as 



Re part: 



"0 



'2 t2,~2 



^Cexp( )cos(— ^ 



aw. 



kH^ 1 

= cosf ) • —= 

+ sm — 

kH'^ 1 

= cosf ) • —=\ 





Its/a jf. I ;iVa ' \ ■"■S 

£/Cexp( jcos — 

awg 2 



^^2^2 2 

li+sin(-— ^)-^Ni 



Because (tildes omitted) 



is satisfied at P, tlie Fresnel approximation is effective at the P (e.g., Ref. 21 §8.3.3). In the present case, setting 
the criteria as <:| ■ ■ ■ | < 2;r x 0.05, the Fresnel approximation is effective in 1 3mm < Z if X < 0.3Z, effective 
in 4.5mm < Z if X < O.IZ, effective in 0.6mm < Z if X = 0. In each case, the moving range of ^ was set to 
-v^vvo < ^ < V2wo. 



44 M.Yasuda 



Impart: / Ji^exp( ^ )sinA;(— 

J-oo Wn 

= —[ 

Vdc J- 



6VS kl 



2s' 



)sm(— T— 



2T2 



272 



= C0S(— r)-^Ni-sin(-— :r)-^Mi . 



where Mi and Ni are 
Mi=/_ 



^-4- ~ (C+-^-jco\/a)^ Ttn 



d^exp(— 



nZ 



)cos— , 



aw; 







)sin — 



^'^p^ H ^ '''' ~ 



Likewise, setting fj = p^its' jk, we write the real and imaginary parts of the fj- 
direction integral as 



Re part: / t/f?exp(- ^^ )cos^ 



(p 



.''II \2 



)cos— = i/— Q, 



Im part: 




in which Q and R represent the integrals of the left sides. For the present purpose, we can 
setyo = 0- Then, because /ewq/tts' is large in the present case, the approximate values of Q 
and R are obtained from the Fresnel integrals: /Q°°cos(7rT^/2)<iT= l/2,/J°sin(a:T^/2)(iT = 
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1 12. That is, Q ~ 1 , R ~ 1 . Thus, Oi is written as (note - #7 4^ 

AS y (X 



) 



U{P) 



jks ' 



^2(1-'-) 

(,c„s(gl, ^ M, +si„(g^) ■ N,) + Kcos(g^) N, - si„(i^) ^ M,) 



-111 



In a same manner, U{P) is obtained as 



2(1-/2) 



/ fe2/'2 j^2J2 j^2J2 jp.J2 

i (cos( — ^) • M2 + sin( — ^) • N2) + /(cos( — ^) • N2 - sin( ^] 



M2) 



where 



M2 = 



''^ _)cos^ 
' 2 



dcwl 



£/^exp(— - 



(C + J^(l-io))^ 



nZ 



kw 



11 
;r2 



-)cos 



2 ' 



+~ . (C- 

^_^flfCexp( 



^-i0\^)2 ^l2 

)sm— ^ 



ii^exp(- 



(C + J^(l-io))^ 



ttZ 



-)sin 



Ml, Ni, M2, and N2 at P = (X,0,Z) on the coordinate system (xjjjz) are expressed 
explicitly, ignoring terms of the second and higher powers of X jZ, as 



M 



exp 



,2 





kw. 



COS c/t 

2 



(I.2a) 



Jl|/'=(X,0,Z) 



exp I 



^0 

;rZ 



sin— — c/t . 
2 



(I.2b) 
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(i+VSP + |z-*))' 



) 



—dx (I.3a) 



exp 



nZ 



COS 



N2|p=(X,0,Z) 




exp 



sin —--dx . 
J 2 



nZ 



(I.3b) 
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